<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "https://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<meta name="generator" content="Doxygen 1.8.16"/>
<meta name="viewport" content="width=device-width, initial-scale=1"/>
<title>Dem Bones: include/DemBones/DemBones.h Source File</title>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="navtree.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="resize.js"></script>
<script type="text/javascript" src="navtreedata.js"></script>
<script type="text/javascript" src="navtree.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
  $(document).ready(initResizable);
/* @license-end */</script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/searchdata.js"></script>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/x-mathjax-config">
  MathJax.Hub.Config({
    extensions: ["tex2jax.js"],
    jax: ["input/TeX","output/HTML-CSS"],
});
</script><script type="text/javascript" async="async" src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js"></script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
 <tbody>
 <tr style="height: 56px;">
  <td id="projectlogo"><img alt="Logo" src="DemBones_small.png"/></td>
  <td id="projectalign" style="padding-left: 0.5em;">
   <div id="projectname">Dem Bones
   &#160;<span id="projectnumber">1.2.0</span>
   </div>
   <div id="projectbrief">Skinning Decomposition Library</div>
  </td>
 </tr>
 </tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.16 -->
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
var searchBox = new SearchBox("searchBox", "search",false,'Search');
/* @license-end */
</script>
<script type="text/javascript" src="menudata.js"></script>
<script type="text/javascript" src="menu.js"></script>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
$(function() {
  initMenu('',true,false,'search.php','Search');
  $(document).ready(function() { init_search(); });
});
/* @license-end */</script>
<div id="main-nav"></div>
</div><!-- top -->
<div id="side-nav" class="ui-resizable side-nav-resizable">
  <div id="nav-tree">
    <div id="nav-tree-contents">
      <div id="nav-sync" class="sync"></div>
    </div>
  </div>
  <div id="splitbar" style="-moz-user-select:none;" 
       class="ui-resizable-handle">
  </div>
</div>
<script type="text/javascript">
/* @license magnet:?xt=urn:btih:cf05388f2679ee054f2beb29a391d25f4e673ac3&amp;dn=gpl-2.0.txt GPL-v2 */
$(document).ready(function(){initNavTree('_dem_bones_8h_source.html','');});
/* @license-end */
</script>
<div id="doc-content">
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
     onmouseover="return searchBox.OnSearchSelectShow()"
     onmouseout="return searchBox.OnSearchSelectHide()"
     onkeydown="return searchBox.OnSearchSelectKey(event)">
</div>

<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0" 
        name="MSearchResults" id="MSearchResults">
</iframe>
</div>

<div class="header">
  <div class="headertitle">
<div class="title">DemBones.h</div>  </div>
</div><!--header-->
<div class="contents">
<a href="_dem_bones_8h.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno">    1</span>&#160;<span class="comment">//               Dem Bones - Skinning Decomposition Library                  //</span></div>
<div class="line"><a name="l00003"></a><span class="lineno">    3</span>&#160;<span class="comment">//         Copyright (c) 2019, Electronic Arts. All rights reserved.         //</span></div>
<div class="line"><a name="l00005"></a><span class="lineno">    5</span>&#160;<span class="comment"></span> </div>
<div class="line"><a name="l00006"></a><span class="lineno">    6</span>&#160; </div>
<div class="line"><a name="l00007"></a><span class="lineno">    7</span>&#160; </div>
<div class="line"><a name="l00008"></a><span class="lineno">    8</span>&#160;<span class="preprocessor">#ifndef DEM_BONES_DEM_BONES</span></div>
<div class="line"><a name="l00009"></a><span class="lineno">    9</span>&#160;<span class="preprocessor">#define DEM_BONES_DEM_BONES</span></div>
<div class="line"><a name="l00010"></a><span class="lineno">   10</span>&#160; </div>
<div class="line"><a name="l00011"></a><span class="lineno">   11</span>&#160;<span class="preprocessor">#include &lt;Eigen/Dense&gt;</span></div>
<div class="line"><a name="l00012"></a><span class="lineno">   12</span>&#160;<span class="preprocessor">#include &lt;Eigen/Sparse&gt;</span></div>
<div class="line"><a name="l00013"></a><span class="lineno">   13</span>&#160;<span class="preprocessor">#include &lt;Eigen/StdVector&gt;</span></div>
<div class="line"><a name="l00014"></a><span class="lineno">   14</span>&#160;<span class="preprocessor">#include &lt;algorithm&gt;</span></div>
<div class="line"><a name="l00015"></a><span class="lineno">   15</span>&#160;<span class="preprocessor">#include &lt;queue&gt;</span></div>
<div class="line"><a name="l00016"></a><span class="lineno">   16</span>&#160;<span class="preprocessor">#include &lt;set&gt;</span></div>
<div class="line"><a name="l00017"></a><span class="lineno">   17</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="_convex_l_s_8h.html">ConvexLS.h</a>&quot;</span></div>
<div class="line"><a name="l00018"></a><span class="lineno">   18</span>&#160; </div>
<div class="line"><a name="l00019"></a><span class="lineno">   19</span>&#160;<span class="preprocessor">#ifndef DEM_BONES_MAT_BLOCKS</span></div>
<div class="line"><a name="l00020"></a><span class="lineno">   20</span>&#160;<span class="preprocessor">#include &quot;<a class="code" href="_mat_blocks_8h.html">MatBlocks.h</a>&quot;</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"><a class="line" href="_dem_bones_8h.html#aaa0daf4013ddd01ce5d2d28e8aaad6f5">   21</a></span>&#160;<span class="preprocessor">#define DEM_BONES_DEM_BONES_MAT_BLOCKS_UNDEFINED</span></div>
<div class="line"><a name="l00022"></a><span class="lineno">   22</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00023"></a><span class="lineno">   23</span>&#160; </div>
<div class="line"><a name="l00024"></a><span class="lineno">   24</span>&#160;<span class="keyword">namespace </span><a class="code" href="namespace_dem.html">Dem</a></div>
<div class="line"><a name="l00025"></a><span class="lineno">   25</span>&#160;{</div>
<div class="line"><a name="l00026"></a><span class="lineno">   26</span>&#160; </div>
<div class="line"><a name="l00065"></a><span class="lineno">   65</span>&#160;<span class="keyword">template</span>&lt;<span class="keyword">class</span> _Scalar, <span class="keyword">class</span> _AniMeshScalar&gt;</div>
<div class="line"><a name="l00066"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html">   66</a></span>&#160;<span class="keyword">class </span><a class="code" href="class_dem_1_1_dem_bones.html">DemBones</a> {</div>
<div class="line"><a name="l00067"></a><span class="lineno">   67</span>&#160;<span class="keyword">public</span>: </div>
<div class="line"><a name="l00068"></a><span class="lineno">   68</span>&#160;    EIGEN_MAKE_ALIGNED_OPERATOR_NEW</div>
<div class="line"><a name="l00069"></a><span class="lineno">   69</span>&#160; </div>
<div class="line"><a name="l00070"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">   70</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a>=Eigen::Matrix&lt;_Scalar, Eigen::Dynamic, Eigen::Dynamic&gt;;</div>
<div class="line"><a name="l00071"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">   71</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Matrix4</a>=Eigen::Matrix&lt;_Scalar, 4, 4&gt;;</div>
<div class="line"><a name="l00072"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aed480cb1bda1b33f20171f7f847c1b06">   72</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#aed480cb1bda1b33f20171f7f847c1b06">Matrix3</a>=Eigen::Matrix&lt;_Scalar, 3, 3&gt;;</div>
<div class="line"><a name="l00073"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">   73</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a>=Eigen::Matrix&lt;_Scalar, Eigen::Dynamic, 1&gt;;</div>
<div class="line"><a name="l00074"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a0c30c98fb363f888345d86e5a11acd14">   74</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#a0c30c98fb363f888345d86e5a11acd14">Vector4</a>=Eigen::Matrix&lt;_Scalar, 4, 1&gt;;</div>
<div class="line"><a name="l00075"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a9404399e2528175da18fd80d7d4fc4f9">   75</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#a9404399e2528175da18fd80d7d4fc4f9">Vector3</a>=Eigen::Matrix&lt;_Scalar, 3, 1&gt;;</div>
<div class="line"><a name="l00076"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aa8ea3a6b166652e954decbf135b3255b">   76</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#aa8ea3a6b166652e954decbf135b3255b">SparseMatrix</a>=Eigen::SparseMatrix&lt;_Scalar&gt;;</div>
<div class="line"><a name="l00077"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">   77</a></span>&#160;    <span class="keyword">using</span> <a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>=Eigen::Triplet&lt;_Scalar&gt;;</div>
<div class="line"><a name="l00078"></a><span class="lineno">   78</span>&#160; </div>
<div class="line"><a name="l00080"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aa7f14b41174dc22e81b089d0c5a56950">   80</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#aa7f14b41174dc22e81b089d0c5a56950">nIters</a>;</div>
<div class="line"><a name="l00081"></a><span class="lineno">   81</span>&#160; </div>
<div class="line"><a name="l00083"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a9079e21634edb0dd2502bfc3a1eedbd0">   83</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#a9079e21634edb0dd2502bfc3a1eedbd0">nInitIters</a>;</div>
<div class="line"><a name="l00084"></a><span class="lineno">   84</span>&#160; </div>
<div class="line"><a name="l00086"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a9afbf2acccc4126e1e76c2a41f3ac649">   86</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#a9afbf2acccc4126e1e76c2a41f3ac649">nTransIters</a>;</div>
<div class="line"><a name="l00088"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a53f9859cd63a8c51f6f50fdfa7251681">   88</a></span>&#160;    _Scalar <a class="code" href="class_dem_1_1_dem_bones.html#a53f9859cd63a8c51f6f50fdfa7251681">transAffine</a>;</div>
<div class="line"><a name="l00090"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a3b125934b440ebcebb3ba44acda41366">   90</a></span>&#160;    _Scalar <a class="code" href="class_dem_1_1_dem_bones.html#a3b125934b440ebcebb3ba44acda41366">transAffineNorm</a>;</div>
<div class="line"><a name="l00091"></a><span class="lineno">   91</span>&#160;    </div>
<div class="line"><a name="l00093"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ae381fbe6610c9f4b169b16106dfbae94">   93</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#ae381fbe6610c9f4b169b16106dfbae94">nWeightsIters</a>;</div>
<div class="line"><a name="l00095"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">   95</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>;</div>
<div class="line"><a name="l00097"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aba10dbe999f33dfb5456c0f3554700a3">   97</a></span>&#160;    _Scalar <a class="code" href="class_dem_1_1_dem_bones.html#aba10dbe999f33dfb5456c0f3554700a3">weightsSmooth</a>;  </div>
<div class="line"><a name="l00099"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a26babdede07e8e220002e420b9161aa9">   99</a></span>&#160;    _Scalar <a class="code" href="class_dem_1_1_dem_bones.html#a26babdede07e8e220002e420b9161aa9">weightsSmoothStep</a>;</div>
<div class="line"><a name="l00101"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">  101</a></span>&#160;    _Scalar <a class="code" href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">weightEps</a>;</div>
<div class="line"><a name="l00102"></a><span class="lineno">  102</span>&#160;    </div>
<div class="line"><a name="l00105"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aa5bb61bfa1a60d355f4cce7f9cc60d5c">  105</a></span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#aa5bb61bfa1a60d355f4cce7f9cc60d5c">DemBones</a>(): <a class="code" href="class_dem_1_1_dem_bones.html#aa7f14b41174dc22e81b089d0c5a56950">nIters</a>(30), <a class="code" href="class_dem_1_1_dem_bones.html#a9079e21634edb0dd2502bfc3a1eedbd0">nInitIters</a>(10),</div>
<div class="line"><a name="l00106"></a><span class="lineno">  106</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a9afbf2acccc4126e1e76c2a41f3ac649">nTransIters</a>(5), <a class="code" href="class_dem_1_1_dem_bones.html#a53f9859cd63a8c51f6f50fdfa7251681">transAffine</a>(_Scalar(10)), <a class="code" href="class_dem_1_1_dem_bones.html#a3b125934b440ebcebb3ba44acda41366">transAffineNorm</a>(_Scalar(4)),</div>
<div class="line"><a name="l00107"></a><span class="lineno">  107</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#ae381fbe6610c9f4b169b16106dfbae94">nWeightsIters</a>(3), <a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>(8), <a class="code" href="class_dem_1_1_dem_bones.html#aba10dbe999f33dfb5456c0f3554700a3">weightsSmooth</a>(_Scalar(1e-4)), <a class="code" href="class_dem_1_1_dem_bones.html#a26babdede07e8e220002e420b9161aa9">weightsSmoothStep</a>(_Scalar(1)),</div>
<div class="line"><a name="l00108"></a><span class="lineno">  108</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">weightEps</a>(_Scalar(1e-15)),</div>
<div class="line"><a name="l00109"></a><span class="lineno">  109</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a06065f12b701c9bb556b7a79d14e984b">iter</a>(_iter), <a class="code" href="class_dem_1_1_dem_bones.html#aed55acae84af67a7e958c830b821e747">iterTransformations</a>(_iterTransformations), <a class="code" href="class_dem_1_1_dem_bones.html#aef1b9e1d0815f32c8f8374376c8d2c22">iterWeights</a>(_iterWeights) {</div>
<div class="line"><a name="l00110"></a><span class="lineno">  110</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a18bcd81e4238a4a298c4ee4ca78ae2fb">clear</a>();</div>
<div class="line"><a name="l00111"></a><span class="lineno">  111</span>&#160;    }</div>
<div class="line"><a name="l00112"></a><span class="lineno">  112</span>&#160;    </div>
<div class="line"><a name="l00114"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">  114</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>;</div>
<div class="line"><a name="l00116"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">  116</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>;                 </div>
<div class="line"><a name="l00118"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">  118</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>;     </div>
<div class="line"><a name="l00120"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">  120</a></span>&#160;    <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; </div>
<div class="line"><a name="l00121"></a><span class="lineno">  121</span>&#160; </div>
<div class="line"><a name="l00123"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#adf3d4f2aa7408bb6a9db540e86fd69c7">  123</a></span>&#160;    Eigen::VectorXi <a class="code" href="class_dem_1_1_dem_bones.html#adf3d4f2aa7408bb6a9db540e86fd69c7">fStart</a>;</div>
<div class="line"><a name="l00125"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">  125</a></span>&#160;    Eigen::VectorXi <a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>;</div>
<div class="line"><a name="l00126"></a><span class="lineno">  126</span>&#160; </div>
<div class="line"><a name="l00128"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">  128</a></span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> <a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>;</div>
<div class="line"><a name="l00129"></a><span class="lineno">  129</span>&#160; </div>
<div class="line"><a name="l00131"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">  131</a></span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#aa8ea3a6b166652e954decbf135b3255b">SparseMatrix</a> <a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>;</div>
<div class="line"><a name="l00132"></a><span class="lineno">  132</span>&#160; </div>
<div class="line"><a name="l00134"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">  134</a></span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> <a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>;</div>
<div class="line"><a name="l00135"></a><span class="lineno">  135</span>&#160; </div>
<div class="line"><a name="l00139"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">  139</a></span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>;</div>
<div class="line"><a name="l00140"></a><span class="lineno">  140</span>&#160; </div>
<div class="line"><a name="l00142"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">  142</a></span>&#160;    Eigen::VectorXi <a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>;</div>
<div class="line"><a name="l00143"></a><span class="lineno">  143</span>&#160; </div>
<div class="line"><a name="l00145"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">  145</a></span>&#160;    Eigen::Matrix&lt;_AniMeshScalar, Eigen::Dynamic, Eigen::Dynamic&gt; <a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>;</div>
<div class="line"><a name="l00146"></a><span class="lineno">  146</span>&#160;    </div>
<div class="line"><a name="l00148"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">  148</a></span>&#160;    std::vector&lt;std::vector&lt;int&gt;&gt; <a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>;</div>
<div class="line"><a name="l00149"></a><span class="lineno">  149</span>&#160; </div>
<div class="line"><a name="l00151"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a06065f12b701c9bb556b7a79d14e984b">  151</a></span>&#160;    <span class="keyword">const</span> <span class="keywordtype">int</span>&amp; <a class="code" href="class_dem_1_1_dem_bones.html#a06065f12b701c9bb556b7a79d14e984b">iter</a>;</div>
<div class="line"><a name="l00152"></a><span class="lineno">  152</span>&#160; </div>
<div class="line"><a name="l00154"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aed55acae84af67a7e958c830b821e747">  154</a></span>&#160;    <span class="keyword">const</span> <span class="keywordtype">int</span>&amp; <a class="code" href="class_dem_1_1_dem_bones.html#aed55acae84af67a7e958c830b821e747">iterTransformations</a>;</div>
<div class="line"><a name="l00155"></a><span class="lineno">  155</span>&#160; </div>
<div class="line"><a name="l00157"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aef1b9e1d0815f32c8f8374376c8d2c22">  157</a></span>&#160;    <span class="keyword">const</span> <span class="keywordtype">int</span>&amp; <a class="code" href="class_dem_1_1_dem_bones.html#aef1b9e1d0815f32c8f8374376c8d2c22">iterWeights</a>;</div>
<div class="line"><a name="l00158"></a><span class="lineno">  158</span>&#160;    </div>
<div class="line"><a name="l00161"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a18bcd81e4238a4a298c4ee4ca78ae2fb">  161</a></span>&#160;    <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a18bcd81e4238a4a298c4ee4ca78ae2fb">clear</a>() {</div>
<div class="line"><a name="l00162"></a><span class="lineno">  162</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>=<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>=<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>=0;</div>
<div class="line"><a name="l00163"></a><span class="lineno">  163</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#adf3d4f2aa7408bb6a9db540e86fd69c7">fStart</a>.resize(0);</div>
<div class="line"><a name="l00164"></a><span class="lineno">  164</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>.resize(0);</div>
<div class="line"><a name="l00165"></a><span class="lineno">  165</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.resize(0, 0);</div>
<div class="line"><a name="l00166"></a><span class="lineno">  166</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.resize(0, 0);</div>
<div class="line"><a name="l00167"></a><span class="lineno">  167</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>.resize(0);</div>
<div class="line"><a name="l00168"></a><span class="lineno">  168</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.resize(0, 0);</div>
<div class="line"><a name="l00169"></a><span class="lineno">  169</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>.resize(0);</div>
<div class="line"><a name="l00170"></a><span class="lineno">  170</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.resize(0, 0);</div>
<div class="line"><a name="l00171"></a><span class="lineno">  171</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>.resize(0);</div>
<div class="line"><a name="l00172"></a><span class="lineno">  172</span>&#160;        modelSize=-1;</div>
<div class="line"><a name="l00173"></a><span class="lineno">  173</span>&#160;        laplacian.resize(0, 0);</div>
<div class="line"><a name="l00174"></a><span class="lineno">  174</span>&#160;    }</div>
<div class="line"><a name="l00175"></a><span class="lineno">  175</span>&#160; </div>
<div class="line"><a name="l00186"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#af1fcf96f35f6c4662fd95998225b9bda">  186</a></span>&#160;    <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#af1fcf96f35f6c4662fd95998225b9bda">init</a>() {</div>
<div class="line"><a name="l00187"></a><span class="lineno">  187</span>&#160;        <span class="keywordflow">if</span> (modelSize&lt;0) modelSize=sqrt((<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>-(<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.rowwise().sum()/<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>).replicate(1, <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>)).squaredNorm()/<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>/<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>);</div>
<div class="line"><a name="l00188"></a><span class="lineno">  188</span>&#160;        <span class="keywordflow">if</span> (laplacian.cols()!=<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>) computeSmoothSolver();</div>
<div class="line"><a name="l00189"></a><span class="lineno">  189</span>&#160; </div>
<div class="line"><a name="l00190"></a><span class="lineno">  190</span>&#160;        <span class="keywordflow">if</span> (((<span class="keywordtype">int</span>)<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.rows()!=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>)||((<span class="keywordtype">int</span>)<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.cols()!=<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>)) { <span class="comment">//No skinning weight</span></div>
<div class="line"><a name="l00191"></a><span class="lineno">  191</span>&#160;            <span class="keywordflow">if</span> (((<span class="keywordtype">int</span>)<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.rows()!=<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>*4)||((<span class="keywordtype">int</span>)<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.cols()!=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4)) { <span class="comment">//No transformation</span></div>
<div class="line"><a name="l00192"></a><span class="lineno">  192</span>&#160;                <span class="keywordtype">int</span> targetNB=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>;</div>
<div class="line"><a name="l00193"></a><span class="lineno">  193</span>&#160;                <span class="comment">//LBG-VQ</span></div>
<div class="line"><a name="l00194"></a><span class="lineno">  194</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>=1;</div>
<div class="line"><a name="l00195"></a><span class="lineno">  195</span>&#160;                label=Eigen::VectorXi::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00196"></a><span class="lineno">  196</span>&#160;                computeTransFromLabel();</div>
<div class="line"><a name="l00197"></a><span class="lineno">  197</span>&#160; </div>
<div class="line"><a name="l00198"></a><span class="lineno">  198</span>&#160;                <span class="keywordtype">bool</span> cont=<span class="keyword">true</span>;</div>
<div class="line"><a name="l00199"></a><span class="lineno">  199</span>&#160;                <span class="keywordflow">while</span> (cont) {</div>
<div class="line"><a name="l00200"></a><span class="lineno">  200</span>&#160;                    <a class="code" href="class_dem_1_1_dem_bones.html#a14f83e468b4a1eb6e2c6c951e7b5a4d3">cbInitSplitBegin</a>();</div>
<div class="line"><a name="l00201"></a><span class="lineno">  201</span>&#160;                    <span class="keywordtype">int</span> prev=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>;</div>
<div class="line"><a name="l00202"></a><span class="lineno">  202</span>&#160;                    split(targetNB, 3);</div>
<div class="line"><a name="l00203"></a><span class="lineno">  203</span>&#160;                    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> rep=0; rep&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a9079e21634edb0dd2502bfc3a1eedbd0">nInitIters</a>; rep++) {</div>
<div class="line"><a name="l00204"></a><span class="lineno">  204</span>&#160;                        computeTransFromLabel();</div>
<div class="line"><a name="l00205"></a><span class="lineno">  205</span>&#160;                        computeLabel();</div>
<div class="line"><a name="l00206"></a><span class="lineno">  206</span>&#160;                        pruneBones(3);</div>
<div class="line"><a name="l00207"></a><span class="lineno">  207</span>&#160;                    }</div>
<div class="line"><a name="l00208"></a><span class="lineno">  208</span>&#160;                    cont=(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>&lt;targetNB)&amp;&amp;(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>&gt;prev);</div>
<div class="line"><a name="l00209"></a><span class="lineno">  209</span>&#160;                    <a class="code" href="class_dem_1_1_dem_bones.html#aeb37e92f6b19ed95f4d2793c5f3dd300">cbInitSplitEnd</a>();</div>
<div class="line"><a name="l00210"></a><span class="lineno">  210</span>&#160;                }</div>
<div class="line"><a name="l00211"></a><span class="lineno">  211</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>=Eigen::VectorXi::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00212"></a><span class="lineno">  212</span>&#160;                labelToWeights();</div>
<div class="line"><a name="l00213"></a><span class="lineno">  213</span>&#160;            } <span class="keywordflow">else</span> initWeights(); <span class="comment">//Has transformations</span></div>
<div class="line"><a name="l00214"></a><span class="lineno">  214</span>&#160;        } <span class="keywordflow">else</span> { <span class="comment">//Has skinning weights</span></div>
<div class="line"><a name="l00215"></a><span class="lineno">  215</span>&#160;            <span class="keywordflow">if</span> (((<span class="keywordtype">int</span>)<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.rows()!=<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>*4)||((<span class="keywordtype">int</span>)<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.cols()!=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4)) { <span class="comment">//No transformation</span></div>
<div class="line"><a name="l00216"></a><span class="lineno">  216</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>=Matrix4::Identity().replicate(<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00217"></a><span class="lineno">  217</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>=Eigen::VectorXi::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00218"></a><span class="lineno">  218</span>&#160;            }</div>
<div class="line"><a name="l00219"></a><span class="lineno">  219</span>&#160;        }</div>
<div class="line"><a name="l00220"></a><span class="lineno">  220</span>&#160; </div>
<div class="line"><a name="l00221"></a><span class="lineno">  221</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>.size()!=<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>) <a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>=VectorX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00222"></a><span class="lineno">  222</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>.size()!=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>) <a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>=Eigen::VectorXi::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00223"></a><span class="lineno">  223</span>&#160;    }</div>
<div class="line"><a name="l00224"></a><span class="lineno">  224</span>&#160; </div>
<div class="line"><a name="l00237"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a673287121735ea64e0e4c244bacc432b">  237</a></span>&#160;    <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a673287121735ea64e0e4c244bacc432b">computeTranformations</a>() {</div>
<div class="line"><a name="l00238"></a><span class="lineno">  238</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#a9afbf2acccc4126e1e76c2a41f3ac649">nTransIters</a>==0) <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00239"></a><span class="lineno">  239</span>&#160; </div>
<div class="line"><a name="l00240"></a><span class="lineno">  240</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#af1fcf96f35f6c4662fd95998225b9bda">init</a>();</div>
<div class="line"><a name="l00241"></a><span class="lineno">  241</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a5d5a50f2d4622a4823065cdca03ac9aa">cbTranformationsBegin</a>();</div>
<div class="line"><a name="l00242"></a><span class="lineno">  242</span>&#160; </div>
<div class="line"><a name="l00243"></a><span class="lineno">  243</span>&#160;        compute_vuT();</div>
<div class="line"><a name="l00244"></a><span class="lineno">  244</span>&#160;        compute_uuT();</div>
<div class="line"><a name="l00245"></a><span class="lineno">  245</span>&#160; </div>
<div class="line"><a name="l00246"></a><span class="lineno">  246</span>&#160;        <span class="keywordflow">for</span> (_iterTransformations=0; _iterTransformations&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a9afbf2acccc4126e1e76c2a41f3ac649">nTransIters</a>; _iterTransformations++) {</div>
<div class="line"><a name="l00247"></a><span class="lineno">  247</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a938532a49b8ae162b866f9b59a75a7ab">cbTransformationsIterBegin</a>();</div>
<div class="line"><a name="l00248"></a><span class="lineno">  248</span>&#160;<span class="preprocessor">            #pragma omp parallel for</span></div>
<div class="line"><a name="l00249"></a><span class="lineno">  249</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++)</div>
<div class="line"><a name="l00250"></a><span class="lineno">  250</span>&#160;                <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) </div>
<div class="line"><a name="l00251"></a><span class="lineno">  251</span>&#160;                    <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">lockM</a>(j)==0) {</div>
<div class="line"><a name="l00252"></a><span class="lineno">  252</span>&#160;                        <a class="code" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Matrix4</a> qpT=vuT.blk4(k, j);</div>
<div class="line"><a name="l00253"></a><span class="lineno">  253</span>&#160;                        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> it=uuT.outerIdx(j); it&lt;uuT.outerIdx(j+1); it++)</div>
<div class="line"><a name="l00254"></a><span class="lineno">  254</span>&#160;                            <span class="keywordflow">if</span> (uuT.innerIdx(it)!=j) qpT-=<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.blk4(k, uuT.innerIdx(it))*uuT.val.blk4(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k), it);</div>
<div class="line"><a name="l00255"></a><span class="lineno">  255</span>&#160;                        qpT2m(qpT, k, j);</div>
<div class="line"><a name="l00256"></a><span class="lineno">  256</span>&#160;                    }</div>
<div class="line"><a name="l00257"></a><span class="lineno">  257</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#aaaf326326d6094637dd0a995a778241f">cbTransformationsIterEnd</a>()) <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00258"></a><span class="lineno">  258</span>&#160;        }</div>
<div class="line"><a name="l00259"></a><span class="lineno">  259</span>&#160;        </div>
<div class="line"><a name="l00260"></a><span class="lineno">  260</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a8008b672232b4d2f4081e00eea2ff062">cbTransformationsEnd</a>();</div>
<div class="line"><a name="l00261"></a><span class="lineno">  261</span>&#160;    }</div>
<div class="line"><a name="l00262"></a><span class="lineno">  262</span>&#160; </div>
<div class="line"><a name="l00276"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a3101275f535d5c2435747b8ea5406f7d">  276</a></span>&#160;    <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a3101275f535d5c2435747b8ea5406f7d">computeWeights</a>() {</div>
<div class="line"><a name="l00277"></a><span class="lineno">  277</span>&#160;        <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#ae381fbe6610c9f4b169b16106dfbae94">nWeightsIters</a>==0) <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00278"></a><span class="lineno">  278</span>&#160;        </div>
<div class="line"><a name="l00279"></a><span class="lineno">  279</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#af1fcf96f35f6c4662fd95998225b9bda">init</a>();</div>
<div class="line"><a name="l00280"></a><span class="lineno">  280</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#ac4698b8e39c770e51aea7122689edcdf">cbWeightsBegin</a>();</div>
<div class="line"><a name="l00281"></a><span class="lineno">  281</span>&#160;        </div>
<div class="line"><a name="l00282"></a><span class="lineno">  282</span>&#160;        compute_mTm();</div>
<div class="line"><a name="l00283"></a><span class="lineno">  283</span>&#160; </div>
<div class="line"><a name="l00284"></a><span class="lineno">  284</span>&#160;        aTb=MatrixX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00285"></a><span class="lineno">  285</span>&#160;        wSolver.init(<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>);</div>
<div class="line"><a name="l00286"></a><span class="lineno">  286</span>&#160;        std::vector&lt;Triplet, Eigen::aligned_allocator&lt;Triplet&gt;&gt; trip;</div>
<div class="line"><a name="l00287"></a><span class="lineno">  287</span>&#160;        trip.reserve(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>*<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>);</div>
<div class="line"><a name="l00288"></a><span class="lineno">  288</span>&#160; </div>
<div class="line"><a name="l00289"></a><span class="lineno">  289</span>&#160;        <span class="keywordflow">for</span> (_iterWeights=0; _iterWeights&lt;<a class="code" href="class_dem_1_1_dem_bones.html#ae381fbe6610c9f4b169b16106dfbae94">nWeightsIters</a>; _iterWeights++) {</div>
<div class="line"><a name="l00290"></a><span class="lineno">  290</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a6391078901046e3acbd582dba6654fdd">cbWeightsIterBegin</a>();</div>
<div class="line"><a name="l00291"></a><span class="lineno">  291</span>&#160; </div>
<div class="line"><a name="l00292"></a><span class="lineno">  292</span>&#160;            compute_ws();</div>
<div class="line"><a name="l00293"></a><span class="lineno">  293</span>&#160;            compute_aTb();</div>
<div class="line"><a name="l00294"></a><span class="lineno">  294</span>&#160; </div>
<div class="line"><a name="l00295"></a><span class="lineno">  295</span>&#160;            <span class="keywordtype">double</span> reg_scale=pow(modelSize, 2)*<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>;</div>
<div class="line"><a name="l00296"></a><span class="lineno">  296</span>&#160; </div>
<div class="line"><a name="l00297"></a><span class="lineno">  297</span>&#160;            trip.clear();</div>
<div class="line"><a name="l00298"></a><span class="lineno">  298</span>&#160;<span class="preprocessor">            #pragma omp parallel for</span></div>
<div class="line"><a name="l00299"></a><span class="lineno">  299</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00300"></a><span class="lineno">  300</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> aTai;</div>
<div class="line"><a name="l00301"></a><span class="lineno">  301</span>&#160;                compute_aTa(i, aTai);</div>
<div class="line"><a name="l00302"></a><span class="lineno">  302</span>&#160;                aTai=(1-<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>(i))*(aTai/reg_scale+<a class="code" href="class_dem_1_1_dem_bones.html#aba10dbe999f33dfb5456c0f3554700a3">weightsSmooth</a>*MatrixX::Identity(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>))+<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>(i)*MatrixX::Identity(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00303"></a><span class="lineno">  303</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> aTbi=(1-<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>(i))*(aTb.col(i)/reg_scale+<a class="code" href="class_dem_1_1_dem_bones.html#aba10dbe999f33dfb5456c0f3554700a3">weightsSmooth</a>*ws.col(i))+<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>(i)*<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.col(i);</div>
<div class="line"><a name="l00304"></a><span class="lineno">  304</span>&#160; </div>
<div class="line"><a name="l00305"></a><span class="lineno">  305</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> x=(1-<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>(i))*ws.col(i)+<a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>(i)*<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.col(i);</div>
<div class="line"><a name="l00306"></a><span class="lineno">  306</span>&#160;                Eigen::ArrayXi idx=Eigen::ArrayXi::LinSpaced(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, 0, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>-1);</div>
<div class="line"><a name="l00307"></a><span class="lineno">  307</span>&#160;                std::sort(idx.data(), idx.data()+<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, [&amp;x](<span class="keywordtype">int</span> i1, <span class="keywordtype">int</span> i2) { <span class="keywordflow">return</span> x(i1)&gt;x(i2); });</div>
<div class="line"><a name="l00308"></a><span class="lineno">  308</span>&#160;                <span class="keywordtype">int</span> nnzi=std::min(<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00309"></a><span class="lineno">  309</span>&#160;                <span class="keywordflow">while</span> (x(idx(nnzi-1))&lt;<a class="code" href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">weightEps</a>) nnzi--;</div>
<div class="line"><a name="l00310"></a><span class="lineno">  310</span>&#160; </div>
<div class="line"><a name="l00311"></a><span class="lineno">  311</span>&#160;                <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> x0=<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.col(i).toDense().cwiseMax(0.0);</div>
<div class="line"><a name="l00312"></a><span class="lineno">  312</span>&#160;                x=<a class="code" href="namespace_dem.html#a5f2dfa58c8c80cb8b494cbdf8d82f105">indexing_vector</a>(x0, idx.head(nnzi));</div>
<div class="line"><a name="l00313"></a><span class="lineno">  313</span>&#160;                _Scalar s=x.sum();</div>
<div class="line"><a name="l00314"></a><span class="lineno">  314</span>&#160;                <span class="keywordflow">if</span> (s&gt;_Scalar(0.1)) x/=s; <span class="keywordflow">else</span> x=VectorX::Constant(nnzi, _Scalar(1)/nnzi);</div>
<div class="line"><a name="l00315"></a><span class="lineno">  315</span>&#160; </div>
<div class="line"><a name="l00316"></a><span class="lineno">  316</span>&#160;                wSolver.solve(<a class="code" href="namespace_dem.html#a86c6e47fcd8031262219c42b65cf84ed">indexing_row_col</a>(aTai, idx.head(nnzi), idx.head(nnzi)), <a class="code" href="namespace_dem.html#a5f2dfa58c8c80cb8b494cbdf8d82f105">indexing_vector</a>(aTbi, idx.head(nnzi)), x, <span class="keyword">true</span>, <span class="keyword">true</span>);</div>
<div class="line"><a name="l00317"></a><span class="lineno">  317</span>&#160; </div>
<div class="line"><a name="l00318"></a><span class="lineno">  318</span>&#160;<span class="preprocessor">                #pragma omp critical</span></div>
<div class="line"><a name="l00319"></a><span class="lineno">  319</span>&#160;                <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;nnzi; j++)</div>
<div class="line"><a name="l00320"></a><span class="lineno">  320</span>&#160;                    <span class="keywordflow">if</span> (x(j)!=0) trip.push_back(<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(idx[j], i, x(j)));</div>
<div class="line"><a name="l00321"></a><span class="lineno">  321</span>&#160;            }</div>
<div class="line"><a name="l00322"></a><span class="lineno">  322</span>&#160; </div>
<div class="line"><a name="l00323"></a><span class="lineno">  323</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.resize(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00324"></a><span class="lineno">  324</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.setFromTriplets(trip.begin(), trip.end());</div>
<div class="line"><a name="l00325"></a><span class="lineno">  325</span>&#160;            </div>
<div class="line"><a name="l00326"></a><span class="lineno">  326</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#a610c7725152df6465cd5f49311a3eeed">cbWeightsIterEnd</a>()) <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00327"></a><span class="lineno">  327</span>&#160;        }</div>
<div class="line"><a name="l00328"></a><span class="lineno">  328</span>&#160;        </div>
<div class="line"><a name="l00329"></a><span class="lineno">  329</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a9dc51d56fe59c8c3e1e221b9819360f9">cbWeightsEnd</a>();</div>
<div class="line"><a name="l00330"></a><span class="lineno">  330</span>&#160;    }</div>
<div class="line"><a name="l00331"></a><span class="lineno">  331</span>&#160; </div>
<div class="line"><a name="l00344"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ac8d2821c730539490ad6ce956a0b62a6">  344</a></span>&#160;    <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#ac8d2821c730539490ad6ce956a0b62a6">compute</a>() {</div>
<div class="line"><a name="l00345"></a><span class="lineno">  345</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#af1fcf96f35f6c4662fd95998225b9bda">init</a>();</div>
<div class="line"><a name="l00346"></a><span class="lineno">  346</span>&#160; </div>
<div class="line"><a name="l00347"></a><span class="lineno">  347</span>&#160;        <span class="keywordflow">for</span> (_iter=0; _iter&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa7f14b41174dc22e81b089d0c5a56950">nIters</a>; _iter++) {</div>
<div class="line"><a name="l00348"></a><span class="lineno">  348</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#af43507c4ae91a4bf6550bc2afc25c023">cbIterBegin</a>();</div>
<div class="line"><a name="l00349"></a><span class="lineno">  349</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a673287121735ea64e0e4c244bacc432b">computeTranformations</a>();</div>
<div class="line"><a name="l00350"></a><span class="lineno">  350</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a3101275f535d5c2435747b8ea5406f7d">computeWeights</a>();</div>
<div class="line"><a name="l00351"></a><span class="lineno">  351</span>&#160;            <span class="keywordflow">if</span> (<a class="code" href="class_dem_1_1_dem_bones.html#a978a27bd67c19d78bd36242a29876212">cbIterEnd</a>()) <span class="keywordflow">break</span>;</div>
<div class="line"><a name="l00352"></a><span class="lineno">  352</span>&#160;        }</div>
<div class="line"><a name="l00353"></a><span class="lineno">  353</span>&#160;    }</div>
<div class="line"><a name="l00354"></a><span class="lineno">  354</span>&#160; </div>
<div class="line"><a name="l00356"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a050b398499f5daa9e42119b0774f25d0">  356</a></span>&#160;    _Scalar <a class="code" href="class_dem_1_1_dem_bones.html#a050b398499f5daa9e42119b0774f25d0">rmse</a>() {</div>
<div class="line"><a name="l00357"></a><span class="lineno">  357</span>&#160;        _Scalar e=0;</div>
<div class="line"><a name="l00358"></a><span class="lineno">  358</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00359"></a><span class="lineno">  359</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00360"></a><span class="lineno">  360</span>&#160;            _Scalar ei=0;</div>
<div class="line"><a name="l00361"></a><span class="lineno">  361</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Matrix4</a> mki;</div>
<div class="line"><a name="l00362"></a><span class="lineno">  362</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++) {</div>
<div class="line"><a name="l00363"></a><span class="lineno">  363</span>&#160;                mki.setZero();</div>
<div class="line"><a name="l00364"></a><span class="lineno">  364</span>&#160;                <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator it(<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>, i); it; ++it) mki+=it.value()*<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.blk4(k, it.row());</div>
<div class="line"><a name="l00365"></a><span class="lineno">  365</span>&#160;                ei+=(mki.template topLeftCorner&lt;3, 3&gt;()*<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k), i)+mki.template topRightCorner&lt;3, 1&gt;()-<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, i).template cast&lt;_Scalar&gt;()).squaredNorm();</div>
<div class="line"><a name="l00366"></a><span class="lineno">  366</span>&#160;            }</div>
<div class="line"><a name="l00367"></a><span class="lineno">  367</span>&#160;<span class="preprocessor">            #pragma omp atomic</span></div>
<div class="line"><a name="l00368"></a><span class="lineno">  368</span>&#160;            e+=ei;</div>
<div class="line"><a name="l00369"></a><span class="lineno">  369</span>&#160;        }</div>
<div class="line"><a name="l00370"></a><span class="lineno">  370</span>&#160;        <span class="keywordflow">return</span> std::sqrt(e/<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>/<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00371"></a><span class="lineno">  371</span>&#160;    }</div>
<div class="line"><a name="l00372"></a><span class="lineno">  372</span>&#160; </div>
<div class="line"><a name="l00374"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a14f83e468b4a1eb6e2c6c951e7b5a4d3">  374</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a14f83e468b4a1eb6e2c6c951e7b5a4d3">cbInitSplitBegin</a>() {}</div>
<div class="line"><a name="l00376"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aeb37e92f6b19ed95f4d2793c5f3dd300">  376</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#aeb37e92f6b19ed95f4d2793c5f3dd300">cbInitSplitEnd</a>() {}</div>
<div class="line"><a name="l00377"></a><span class="lineno">  377</span>&#160; </div>
<div class="line"><a name="l00379"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#af43507c4ae91a4bf6550bc2afc25c023">  379</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#af43507c4ae91a4bf6550bc2afc25c023">cbIterBegin</a>() {}</div>
<div class="line"><a name="l00381"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a978a27bd67c19d78bd36242a29876212">  381</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="class_dem_1_1_dem_bones.html#a978a27bd67c19d78bd36242a29876212">cbIterEnd</a>() { <span class="keywordflow">return</span> <span class="keyword">false</span>; }</div>
<div class="line"><a name="l00382"></a><span class="lineno">  382</span>&#160; </div>
<div class="line"><a name="l00384"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#ac4698b8e39c770e51aea7122689edcdf">  384</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#ac4698b8e39c770e51aea7122689edcdf">cbWeightsBegin</a>() {}</div>
<div class="line"><a name="l00386"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a9dc51d56fe59c8c3e1e221b9819360f9">  386</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a9dc51d56fe59c8c3e1e221b9819360f9">cbWeightsEnd</a>() {}</div>
<div class="line"><a name="l00387"></a><span class="lineno">  387</span>&#160; </div>
<div class="line"><a name="l00389"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a5d5a50f2d4622a4823065cdca03ac9aa">  389</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a5d5a50f2d4622a4823065cdca03ac9aa">cbTranformationsBegin</a>() {}</div>
<div class="line"><a name="l00391"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a8008b672232b4d2f4081e00eea2ff062">  391</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a8008b672232b4d2f4081e00eea2ff062">cbTransformationsEnd</a>() {}</div>
<div class="line"><a name="l00392"></a><span class="lineno">  392</span>&#160; </div>
<div class="line"><a name="l00394"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a938532a49b8ae162b866f9b59a75a7ab">  394</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a938532a49b8ae162b866f9b59a75a7ab">cbTransformationsIterBegin</a>() {}</div>
<div class="line"><a name="l00396"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#aaaf326326d6094637dd0a995a778241f">  396</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="class_dem_1_1_dem_bones.html#aaaf326326d6094637dd0a995a778241f">cbTransformationsIterEnd</a>() { <span class="keywordflow">return</span> <span class="keyword">false</span>; }</div>
<div class="line"><a name="l00397"></a><span class="lineno">  397</span>&#160; </div>
<div class="line"><a name="l00399"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a6391078901046e3acbd582dba6654fdd">  399</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">void</span> <a class="code" href="class_dem_1_1_dem_bones.html#a6391078901046e3acbd582dba6654fdd">cbWeightsIterBegin</a>() {}</div>
<div class="line"><a name="l00401"></a><span class="lineno"><a class="line" href="class_dem_1_1_dem_bones.html#a610c7725152df6465cd5f49311a3eeed">  401</a></span>&#160;    <span class="keyword">virtual</span> <span class="keywordtype">bool</span> <a class="code" href="class_dem_1_1_dem_bones.html#a610c7725152df6465cd5f49311a3eeed">cbWeightsIterEnd</a>() { <span class="keywordflow">return</span> <span class="keyword">false</span>; }</div>
<div class="line"><a name="l00402"></a><span class="lineno">  402</span>&#160; </div>
<div class="line"><a name="l00403"></a><span class="lineno">  403</span>&#160;<span class="keyword">private</span>:</div>
<div class="line"><a name="l00404"></a><span class="lineno">  404</span>&#160;    <span class="keywordtype">int</span> _iter, _iterTransformations, _iterWeights;</div>
<div class="line"><a name="l00405"></a><span class="lineno">  405</span>&#160; </div>
<div class="line"><a name="l00411"></a><span class="lineno">  411</span>&#160;    <span class="keywordtype">void</span> qpT2m(<span class="keyword">const</span> <a class="code" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Matrix4</a>&amp; _qpT, <span class="keywordtype">int</span> k, <span class="keywordtype">int</span> j) {</div>
<div class="line"><a name="l00412"></a><span class="lineno">  412</span>&#160;        <span class="keywordflow">if</span> (_qpT(3, 3)!=0) {</div>
<div class="line"><a name="l00413"></a><span class="lineno">  413</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Matrix4</a> qpT=_qpT/_qpT(3, 3);</div>
<div class="line"><a name="l00414"></a><span class="lineno">  414</span>&#160;            Eigen::JacobiSVD&lt;Matrix3&gt; svd(qpT.template topLeftCorner&lt;3, 3&gt;()-qpT.template topRightCorner&lt;3, 1&gt;()*qpT.template bottomLeftCorner&lt;1, 3&gt;(), Eigen::ComputeFullU|Eigen::ComputeFullV);</div>
<div class="line"><a name="l00415"></a><span class="lineno">  415</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#aed480cb1bda1b33f20171f7f847c1b06">Matrix3</a> d=Matrix3::Identity();</div>
<div class="line"><a name="l00416"></a><span class="lineno">  416</span>&#160;            d(2, 2)=(svd.matrixU()*svd.matrixV().transpose()).determinant();</div>
<div class="line"><a name="l00417"></a><span class="lineno">  417</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.rotMat(k, j)=svd.matrixU()*d*svd.matrixV().transpose();</div>
<div class="line"><a name="l00418"></a><span class="lineno">  418</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.transVec(k, j)=qpT.template topRightCorner&lt;3, 1&gt;()-<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.rotMat(k, j)*qpT.template bottomLeftCorner&lt;1, 3&gt;().transpose();</div>
<div class="line"><a name="l00419"></a><span class="lineno">  419</span>&#160;        }</div>
<div class="line"><a name="l00420"></a><span class="lineno">  420</span>&#160;    }</div>
<div class="line"><a name="l00421"></a><span class="lineno">  421</span>&#160; </div>
<div class="line"><a name="l00426"></a><span class="lineno">  426</span>&#160;    _Scalar errorVtxBone(<span class="keywordtype">int</span> i, <span class="keywordtype">int</span> j, <span class="keywordtype">bool</span> par=<span class="keyword">true</span>) {</div>
<div class="line"><a name="l00427"></a><span class="lineno">  427</span>&#160;        _Scalar e=0;</div>
<div class="line"><a name="l00428"></a><span class="lineno">  428</span>&#160;<span class="preprocessor">        #pragma omp parallel for if(par)</span></div>
<div class="line"><a name="l00429"></a><span class="lineno">  429</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++)</div>
<div class="line"><a name="l00430"></a><span class="lineno">  430</span>&#160;            #pragma omp atomic</div>
<div class="line"><a name="l00431"></a><span class="lineno">  431</span>&#160;            e+=(<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.rotMat(k, j)*<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k), i)+<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.transVec(k, j)-<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, i).template cast&lt;_Scalar&gt;()).squaredNorm();</div>
<div class="line"><a name="l00432"></a><span class="lineno">  432</span>&#160;        <span class="keywordflow">return</span> e;</div>
<div class="line"><a name="l00433"></a><span class="lineno">  433</span>&#160;    }</div>
<div class="line"><a name="l00434"></a><span class="lineno">  434</span>&#160; </div>
<div class="line"><a name="l00436"></a><span class="lineno">  436</span>&#160;    Eigen::VectorXi label;</div>
<div class="line"><a name="l00437"></a><span class="lineno">  437</span>&#160; </div>
<div class="line"><a name="l00439"></a><span class="lineno">  439</span>&#160;    <span class="keyword">struct </span>TripletLess {</div>
<div class="line"><a name="l00440"></a><span class="lineno">  440</span>&#160;        <span class="keywordtype">bool</span> operator() (<span class="keyword">const</span> <a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>&amp; t1, <span class="keyword">const</span> <a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>&amp; t2) {</div>
<div class="line"><a name="l00441"></a><span class="lineno">  441</span>&#160;            <span class="keywordflow">return</span> t1.value()&gt;t2.value();</div>
<div class="line"><a name="l00442"></a><span class="lineno">  442</span>&#160;        }</div>
<div class="line"><a name="l00443"></a><span class="lineno">  443</span>&#160;    };</div>
<div class="line"><a name="l00444"></a><span class="lineno">  444</span>&#160; </div>
<div class="line"><a name="l00447"></a><span class="lineno">  447</span>&#160;    <span class="keywordtype">void</span> computeLabel() {</div>
<div class="line"><a name="l00448"></a><span class="lineno">  448</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> ei(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00449"></a><span class="lineno">  449</span>&#160;        Eigen::VectorXi seed=Eigen::VectorXi::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, -1);</div>
<div class="line"><a name="l00450"></a><span class="lineno">  450</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> gMin(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00451"></a><span class="lineno">  451</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00452"></a><span class="lineno">  452</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00453"></a><span class="lineno">  453</span>&#160;            <span class="keywordtype">int</span> j=label(i);</div>
<div class="line"><a name="l00454"></a><span class="lineno">  454</span>&#160;            <span class="keywordflow">if</span> (j!=-1) {</div>
<div class="line"><a name="l00455"></a><span class="lineno">  455</span>&#160;                ei(i)=errorVtxBone(i, j, <span class="keyword">false</span>);</div>
<div class="line"><a name="l00456"></a><span class="lineno">  456</span>&#160;                <span class="keywordflow">if</span> ((seed(j)==-1)||(ei(i)&lt;gMin(j))) {</div>
<div class="line"><a name="l00457"></a><span class="lineno">  457</span>&#160;<span class="preprocessor">                    #pragma omp critical</span></div>
<div class="line"><a name="l00458"></a><span class="lineno">  458</span>&#160;                    <span class="keywordflow">if</span> ((seed(j)==-1)||(ei(i)&lt;gMin(j))) {</div>
<div class="line"><a name="l00459"></a><span class="lineno">  459</span>&#160;                        gMin(j)=ei(i);</div>
<div class="line"><a name="l00460"></a><span class="lineno">  460</span>&#160;                        seed(j)=i;</div>
<div class="line"><a name="l00461"></a><span class="lineno">  461</span>&#160;                    }</div>
<div class="line"><a name="l00462"></a><span class="lineno">  462</span>&#160;                }</div>
<div class="line"><a name="l00463"></a><span class="lineno">  463</span>&#160;            }</div>
<div class="line"><a name="l00464"></a><span class="lineno">  464</span>&#160;        }</div>
<div class="line"><a name="l00465"></a><span class="lineno">  465</span>&#160; </div>
<div class="line"><a name="l00466"></a><span class="lineno">  466</span>&#160;        std::priority_queue&lt;Triplet, std::vector&lt;Triplet, Eigen::aligned_allocator&lt;Triplet&gt;&gt;, TripletLess&gt; heap;</div>
<div class="line"><a name="l00467"></a><span class="lineno">  467</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) <span class="keywordflow">if</span> (seed(j)!=-1) heap.push(<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(j, seed(j), ei(seed(j))));</div>
<div class="line"><a name="l00468"></a><span class="lineno">  468</span>&#160; </div>
<div class="line"><a name="l00469"></a><span class="lineno">  469</span>&#160;        <span class="keywordflow">if</span> (laplacian.cols()!=<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>) computeSmoothSolver();</div>
<div class="line"><a name="l00470"></a><span class="lineno">  470</span>&#160; </div>
<div class="line"><a name="l00471"></a><span class="lineno">  471</span>&#160;        std::vector&lt;bool&gt; dirty(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>, <span class="keyword">true</span>);</div>
<div class="line"><a name="l00472"></a><span class="lineno">  472</span>&#160;        <span class="keywordflow">while</span> (!heap.empty()) {</div>
<div class="line"><a name="l00473"></a><span class="lineno">  473</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a> top=heap.top();</div>
<div class="line"><a name="l00474"></a><span class="lineno">  474</span>&#160;            heap.pop();</div>
<div class="line"><a name="l00475"></a><span class="lineno">  475</span>&#160;            <span class="keywordtype">int</span> i=(int)top.col();</div>
<div class="line"><a name="l00476"></a><span class="lineno">  476</span>&#160;            <span class="keywordtype">int</span> j=(int)top.row();</div>
<div class="line"><a name="l00477"></a><span class="lineno">  477</span>&#160;            <span class="keywordflow">if</span> (dirty[i]) {</div>
<div class="line"><a name="l00478"></a><span class="lineno">  478</span>&#160;                label(i)=j;</div>
<div class="line"><a name="l00479"></a><span class="lineno">  479</span>&#160;                ei(i)=top.value();</div>
<div class="line"><a name="l00480"></a><span class="lineno">  480</span>&#160;                dirty[i]=<span class="keyword">false</span>;</div>
<div class="line"><a name="l00481"></a><span class="lineno">  481</span>&#160;                <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator it(laplacian, i); it; ++it) {</div>
<div class="line"><a name="l00482"></a><span class="lineno">  482</span>&#160;                    <span class="keywordtype">int</span> i2=(int)it.row();</div>
<div class="line"><a name="l00483"></a><span class="lineno">  483</span>&#160;                    <span class="keywordflow">if</span> (dirty[i2]) {</div>
<div class="line"><a name="l00484"></a><span class="lineno">  484</span>&#160;                        <span class="keywordtype">double</span> tmp=(label(i2)==j)?ei(i2):errorVtxBone(i2, j);</div>
<div class="line"><a name="l00485"></a><span class="lineno">  485</span>&#160;                        heap.push(<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(j, i2, tmp));</div>
<div class="line"><a name="l00486"></a><span class="lineno">  486</span>&#160;                    }</div>
<div class="line"><a name="l00487"></a><span class="lineno">  487</span>&#160;                }</div>
<div class="line"><a name="l00488"></a><span class="lineno">  488</span>&#160;            }</div>
<div class="line"><a name="l00489"></a><span class="lineno">  489</span>&#160;        }</div>
<div class="line"><a name="l00490"></a><span class="lineno">  490</span>&#160; </div>
<div class="line"><a name="l00491"></a><span class="lineno">  491</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00492"></a><span class="lineno">  492</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) </div>
<div class="line"><a name="l00493"></a><span class="lineno">  493</span>&#160;            <span class="keywordflow">if</span> (label(i)==-1) {</div>
<div class="line"><a name="l00494"></a><span class="lineno">  494</span>&#160;                _Scalar gMin;</div>
<div class="line"><a name="l00495"></a><span class="lineno">  495</span>&#160;                <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) {</div>
<div class="line"><a name="l00496"></a><span class="lineno">  496</span>&#160;                    _Scalar ej=errorVtxBone(i, j, <span class="keyword">false</span>);</div>
<div class="line"><a name="l00497"></a><span class="lineno">  497</span>&#160;                    <span class="keywordflow">if</span> ((label(i)==-1)||(gMin&gt;ej)) {</div>
<div class="line"><a name="l00498"></a><span class="lineno">  498</span>&#160;                        gMin=ej;</div>
<div class="line"><a name="l00499"></a><span class="lineno">  499</span>&#160;                        label(i)=j;</div>
<div class="line"><a name="l00500"></a><span class="lineno">  500</span>&#160;                    }</div>
<div class="line"><a name="l00501"></a><span class="lineno">  501</span>&#160;                }</div>
<div class="line"><a name="l00502"></a><span class="lineno">  502</span>&#160;            }</div>
<div class="line"><a name="l00503"></a><span class="lineno">  503</span>&#160;    }</div>
<div class="line"><a name="l00504"></a><span class="lineno">  504</span>&#160; </div>
<div class="line"><a name="l00507"></a><span class="lineno">  507</span>&#160;    <span class="keywordtype">void</span> computeTransFromLabel() {</div>
<div class="line"><a name="l00508"></a><span class="lineno">  508</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>=Matrix4::Identity().replicate(<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00509"></a><span class="lineno">  509</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00510"></a><span class="lineno">  510</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++) {</div>
<div class="line"><a name="l00511"></a><span class="lineno">  511</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> qpT=MatrixX::Zero(4, 4*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00512"></a><span class="lineno">  512</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) </div>
<div class="line"><a name="l00513"></a><span class="lineno">  513</span>&#160;                <span class="keywordflow">if</span> (label(i)!=-1) qpT.blk4(0, label(i))+=<a class="code" href="class_dem_1_1_dem_bones.html#a0c30c98fb363f888345d86e5a11acd14">Vector4</a>(<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, i).template cast&lt;_Scalar&gt;().homogeneous())*<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k), i).homogeneous().transpose();</div>
<div class="line"><a name="l00514"></a><span class="lineno">  514</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) qpT2m(qpT.blk4(0, j), k, j);</div>
<div class="line"><a name="l00515"></a><span class="lineno">  515</span>&#160;        }</div>
<div class="line"><a name="l00516"></a><span class="lineno">  516</span>&#160;    }</div>
<div class="line"><a name="l00517"></a><span class="lineno">  517</span>&#160; </div>
<div class="line"><a name="l00520"></a><span class="lineno">  520</span>&#160;    <span class="keywordtype">void</span> labelToWeights() {</div>
<div class="line"><a name="l00521"></a><span class="lineno">  521</span>&#160;        std::vector&lt;Triplet, Eigen::aligned_allocator&lt;Triplet&gt;&gt; trip(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00522"></a><span class="lineno">  522</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) trip[i]=<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(label(i), i, _Scalar(1));</div>
<div class="line"><a name="l00523"></a><span class="lineno">  523</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.resize(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00524"></a><span class="lineno">  524</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.setFromTriplets(trip.begin(), trip.end());</div>
<div class="line"><a name="l00525"></a><span class="lineno">  525</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">lockW</a>=VectorX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00526"></a><span class="lineno">  526</span>&#160;    }</div>
<div class="line"><a name="l00527"></a><span class="lineno">  527</span>&#160; </div>
<div class="line"><a name="l00532"></a><span class="lineno">  532</span>&#160;    <span class="keywordtype">void</span> split(<span class="keywordtype">int</span> maxB, <span class="keywordtype">int</span> threshold) {</div>
<div class="line"><a name="l00533"></a><span class="lineno">  533</span>&#160;        <span class="comment">//Centroids</span></div>
<div class="line"><a name="l00534"></a><span class="lineno">  534</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> cu=MatrixX::Zero(3*<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00535"></a><span class="lineno">  535</span>&#160;        Eigen::VectorXi s=Eigen::VectorXi::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00536"></a><span class="lineno">  536</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00537"></a><span class="lineno">  537</span>&#160;            cu.col(label(i))+=<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.col(i);</div>
<div class="line"><a name="l00538"></a><span class="lineno">  538</span>&#160;            s(label(i))++;</div>
<div class="line"><a name="l00539"></a><span class="lineno">  539</span>&#160;        }</div>
<div class="line"><a name="l00540"></a><span class="lineno">  540</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) <span class="keywordflow">if</span> (s(j)!=0) cu.col(j)/=_Scalar(s(j));</div>
<div class="line"><a name="l00541"></a><span class="lineno">  541</span>&#160; </div>
<div class="line"><a name="l00542"></a><span class="lineno">  542</span>&#160;        <span class="comment">//Distance to centroid &amp; error</span></div>
<div class="line"><a name="l00543"></a><span class="lineno">  543</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> d(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>), e(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00544"></a><span class="lineno">  544</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> minD=VectorX::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, std::numeric_limits&lt;_Scalar&gt;::max());</div>
<div class="line"><a name="l00545"></a><span class="lineno">  545</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> minE=VectorX::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, std::numeric_limits&lt;_Scalar&gt;::max());</div>
<div class="line"><a name="l00546"></a><span class="lineno">  546</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> ce=VectorX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00547"></a><span class="lineno">  547</span>&#160; </div>
<div class="line"><a name="l00548"></a><span class="lineno">  548</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00549"></a><span class="lineno">  549</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00550"></a><span class="lineno">  550</span>&#160;            <span class="keywordtype">int</span> j=label(i);</div>
<div class="line"><a name="l00551"></a><span class="lineno">  551</span>&#160;            d(i)=(<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.col(i)-cu.col(j)).norm();</div>
<div class="line"><a name="l00552"></a><span class="lineno">  552</span>&#160;            e(i)=sqrt(errorVtxBone(i, j, <span class="keyword">false</span>));</div>
<div class="line"><a name="l00553"></a><span class="lineno">  553</span>&#160;            <span class="keywordflow">if</span> (d(i)&lt;minD(j)) {</div>
<div class="line"><a name="l00554"></a><span class="lineno">  554</span>&#160;<span class="preprocessor">                #pragma omp critical</span></div>
<div class="line"><a name="l00555"></a><span class="lineno">  555</span>&#160;                minD(j)=std::min(minD(j), d(i));</div>
<div class="line"><a name="l00556"></a><span class="lineno">  556</span>&#160;            }</div>
<div class="line"><a name="l00557"></a><span class="lineno">  557</span>&#160;            <span class="keywordflow">if</span> (e(i)&lt;minE(j)) {</div>
<div class="line"><a name="l00558"></a><span class="lineno">  558</span>&#160;<span class="preprocessor">                #pragma omp critical</span></div>
<div class="line"><a name="l00559"></a><span class="lineno">  559</span>&#160;                minE(j)=std::min(minE(j), e(i));</div>
<div class="line"><a name="l00560"></a><span class="lineno">  560</span>&#160;            }</div>
<div class="line"><a name="l00561"></a><span class="lineno">  561</span>&#160;<span class="preprocessor">            #pragma omp atomic</span></div>
<div class="line"><a name="l00562"></a><span class="lineno">  562</span>&#160;            ce(j)+=e(i);</div>
<div class="line"><a name="l00563"></a><span class="lineno">  563</span>&#160;        }</div>
<div class="line"><a name="l00564"></a><span class="lineno">  564</span>&#160; </div>
<div class="line"><a name="l00565"></a><span class="lineno">  565</span>&#160;        <span class="comment">//Seed</span></div>
<div class="line"><a name="l00566"></a><span class="lineno">  566</span>&#160;        Eigen::VectorXi seed=Eigen::VectorXi::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, -1);</div>
<div class="line"><a name="l00567"></a><span class="lineno">  567</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> gMax(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00568"></a><span class="lineno">  568</span>&#160; </div>
<div class="line"><a name="l00569"></a><span class="lineno">  569</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00570"></a><span class="lineno">  570</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00571"></a><span class="lineno">  571</span>&#160;            <span class="keywordtype">int</span> j=label(i);</div>
<div class="line"><a name="l00572"></a><span class="lineno">  572</span>&#160;            <span class="keywordtype">double</span> tmp=abs((e(i)-minE(j))*(d(i)-minD(j)));</div>
<div class="line"><a name="l00573"></a><span class="lineno">  573</span>&#160; </div>
<div class="line"><a name="l00574"></a><span class="lineno">  574</span>&#160;            <span class="keywordflow">if</span> ((seed(j)==-1)||(tmp&gt;gMax(j))) {</div>
<div class="line"><a name="l00575"></a><span class="lineno">  575</span>&#160;<span class="preprocessor">                #pragma omp critical</span></div>
<div class="line"><a name="l00576"></a><span class="lineno">  576</span>&#160;                <span class="keywordflow">if</span> ((seed(j)==-1)||(tmp&gt;gMax(j))) {</div>
<div class="line"><a name="l00577"></a><span class="lineno">  577</span>&#160;                    gMax(j)=tmp;</div>
<div class="line"><a name="l00578"></a><span class="lineno">  578</span>&#160;                    seed(j)=i;</div>
<div class="line"><a name="l00579"></a><span class="lineno">  579</span>&#160;                }</div>
<div class="line"><a name="l00580"></a><span class="lineno">  580</span>&#160;            }</div>
<div class="line"><a name="l00581"></a><span class="lineno">  581</span>&#160;        }</div>
<div class="line"><a name="l00582"></a><span class="lineno">  582</span>&#160; </div>
<div class="line"><a name="l00583"></a><span class="lineno">  583</span>&#160;        <span class="keywordtype">int</span> countID=<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>;</div>
<div class="line"><a name="l00584"></a><span class="lineno">  584</span>&#160;        _Scalar avgErr=ce.sum()/<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>;</div>
<div class="line"><a name="l00585"></a><span class="lineno">  585</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++)</div>
<div class="line"><a name="l00586"></a><span class="lineno">  586</span>&#160;            <span class="keywordflow">if</span> ((countID&lt;maxB)&amp;&amp;(s(j)&gt;threshold*2)&amp;&amp;(ce(j)&gt;avgErr/100)) {</div>
<div class="line"><a name="l00587"></a><span class="lineno">  587</span>&#160;                <span class="keywordtype">int</span> newLabel=countID++;</div>
<div class="line"><a name="l00588"></a><span class="lineno">  588</span>&#160;                <span class="keywordtype">int</span> i=seed(j);</div>
<div class="line"><a name="l00589"></a><span class="lineno">  589</span>&#160;                <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator it(laplacian, i); it; ++it) label(it.row())=newLabel;</div>
<div class="line"><a name="l00590"></a><span class="lineno">  590</span>&#160;            }</div>
<div class="line"><a name="l00591"></a><span class="lineno">  591</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>=countID;</div>
<div class="line"><a name="l00592"></a><span class="lineno">  592</span>&#160;    }</div>
<div class="line"><a name="l00593"></a><span class="lineno">  593</span>&#160; </div>
<div class="line"><a name="l00597"></a><span class="lineno">  597</span>&#160;    <span class="keywordtype">void</span> pruneBones(<span class="keywordtype">int</span> threshold) {</div>
<div class="line"><a name="l00598"></a><span class="lineno">  598</span>&#160;        Eigen::VectorXi s=Eigen::VectorXi::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00599"></a><span class="lineno">  599</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00600"></a><span class="lineno">  600</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00601"></a><span class="lineno">  601</span>&#160;<span class="preprocessor">            #pragma omp atomic</span></div>
<div class="line"><a name="l00602"></a><span class="lineno">  602</span>&#160;            s(label(i))++;</div>
<div class="line"><a name="l00603"></a><span class="lineno">  603</span>&#160;        }</div>
<div class="line"><a name="l00604"></a><span class="lineno">  604</span>&#160; </div>
<div class="line"><a name="l00605"></a><span class="lineno">  605</span>&#160;        Eigen::VectorXi newID(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00606"></a><span class="lineno">  606</span>&#160;        <span class="keywordtype">int</span> countID=0;</div>
<div class="line"><a name="l00607"></a><span class="lineno">  607</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++)</div>
<div class="line"><a name="l00608"></a><span class="lineno">  608</span>&#160;            <span class="keywordflow">if</span> (s(j)&lt;threshold) newID(j)=-1; <span class="keywordflow">else</span> newID(j)=countID++;</div>
<div class="line"><a name="l00609"></a><span class="lineno">  609</span>&#160; </div>
<div class="line"><a name="l00610"></a><span class="lineno">  610</span>&#160;        <span class="keywordflow">if</span> (countID==<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>) <span class="keywordflow">return</span>;</div>
<div class="line"><a name="l00611"></a><span class="lineno">  611</span>&#160; </div>
<div class="line"><a name="l00612"></a><span class="lineno">  612</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++)</div>
<div class="line"><a name="l00613"></a><span class="lineno">  613</span>&#160;            <span class="keywordflow">if</span> (newID(j)!=-1) <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.template middleCols&lt;4&gt;(newID(j)*4)=<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.template middleCols&lt;4&gt;(j*4);</div>
<div class="line"><a name="l00614"></a><span class="lineno">  614</span>&#160; </div>
<div class="line"><a name="l00615"></a><span class="lineno">  615</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00616"></a><span class="lineno">  616</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) label(i)=newID(label(i));</div>
<div class="line"><a name="l00617"></a><span class="lineno">  617</span>&#160; </div>
<div class="line"><a name="l00618"></a><span class="lineno">  618</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>=countID;</div>
<div class="line"><a name="l00619"></a><span class="lineno">  619</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.conservativeResize(<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>*4, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4);</div>
<div class="line"><a name="l00620"></a><span class="lineno">  620</span>&#160;        computeLabel();</div>
<div class="line"><a name="l00621"></a><span class="lineno">  621</span>&#160;    }</div>
<div class="line"><a name="l00622"></a><span class="lineno">  622</span>&#160; </div>
<div class="line"><a name="l00625"></a><span class="lineno">  625</span>&#160;    <span class="keywordtype">void</span> initWeights() {</div>
<div class="line"><a name="l00626"></a><span class="lineno">  626</span>&#160;        label=Eigen::VectorXi::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>, -1);</div>
<div class="line"><a name="l00627"></a><span class="lineno">  627</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00628"></a><span class="lineno">  628</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00629"></a><span class="lineno">  629</span>&#160;            _Scalar gMin;</div>
<div class="line"><a name="l00630"></a><span class="lineno">  630</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) {</div>
<div class="line"><a name="l00631"></a><span class="lineno">  631</span>&#160;                _Scalar ej=errorVtxBone(i, j, <span class="keyword">false</span>);</div>
<div class="line"><a name="l00632"></a><span class="lineno">  632</span>&#160;                <span class="keywordflow">if</span> ((label(i)==-1)||(gMin&gt;ej)) {</div>
<div class="line"><a name="l00633"></a><span class="lineno">  633</span>&#160;                    gMin=ej;</div>
<div class="line"><a name="l00634"></a><span class="lineno">  634</span>&#160;                    label(i)=j;</div>
<div class="line"><a name="l00635"></a><span class="lineno">  635</span>&#160;                }</div>
<div class="line"><a name="l00636"></a><span class="lineno">  636</span>&#160;            }</div>
<div class="line"><a name="l00637"></a><span class="lineno">  637</span>&#160;        }</div>
<div class="line"><a name="l00638"></a><span class="lineno">  638</span>&#160;        computeLabel();</div>
<div class="line"><a name="l00639"></a><span class="lineno">  639</span>&#160;        labelToWeights();</div>
<div class="line"><a name="l00640"></a><span class="lineno">  640</span>&#160;    }</div>
<div class="line"><a name="l00641"></a><span class="lineno">  641</span>&#160; </div>
<div class="line"><a name="l00643"></a><span class="lineno">  643</span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> vuT;</div>
<div class="line"><a name="l00644"></a><span class="lineno">  644</span>&#160; </div>
<div class="line"><a name="l00647"></a><span class="lineno">  647</span>&#160;    <span class="keywordtype">void</span> compute_vuT() {</div>
<div class="line"><a name="l00648"></a><span class="lineno">  648</span>&#160;        vuT=MatrixX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>*4, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4);</div>
<div class="line"><a name="l00649"></a><span class="lineno">  649</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00650"></a><span class="lineno">  650</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++) {</div>
<div class="line"><a name="l00651"></a><span class="lineno">  651</span>&#160;            <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> vuTp=MatrixX::Zero(4, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4);</div>
<div class="line"><a name="l00652"></a><span class="lineno">  652</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++)</div>
<div class="line"><a name="l00653"></a><span class="lineno">  653</span>&#160;                <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator it(<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>, i); it; ++it) {</div>
<div class="line"><a name="l00654"></a><span class="lineno">  654</span>&#160;                    <a class="code" href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Matrix4</a> tmp=<a class="code" href="class_dem_1_1_dem_bones.html#a0c30c98fb363f888345d86e5a11acd14">Vector4</a>(<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, i).template cast&lt;_Scalar&gt;().homogeneous())*<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k), i).homogeneous().transpose();</div>
<div class="line"><a name="l00655"></a><span class="lineno">  655</span>&#160;                    vuT.blk4(k, it.row())+=it.value()*tmp;</div>
<div class="line"><a name="l00656"></a><span class="lineno">  656</span>&#160;                    vuTp.blk4(0, it.row())+=pow(it.value(), <a class="code" href="class_dem_1_1_dem_bones.html#a3b125934b440ebcebb3ba44acda41366">transAffineNorm</a>)*tmp;</div>
<div class="line"><a name="l00657"></a><span class="lineno">  657</span>&#160;                }</div>
<div class="line"><a name="l00658"></a><span class="lineno">  658</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++)</div>
<div class="line"><a name="l00659"></a><span class="lineno">  659</span>&#160;                <span class="keywordflow">if</span> (vuTp(3, j*4+3)!=0)</div>
<div class="line"><a name="l00660"></a><span class="lineno">  660</span>&#160;                    vuT.blk4(k, j)+=(<a class="code" href="class_dem_1_1_dem_bones.html#a53f9859cd63a8c51f6f50fdfa7251681">transAffine</a>*vuT(k*4+3, j*4+3)/vuTp(3, j*4+3))*vuTp.blk4(0, j);</div>
<div class="line"><a name="l00661"></a><span class="lineno">  661</span>&#160;        }</div>
<div class="line"><a name="l00662"></a><span class="lineno">  662</span>&#160;    }</div>
<div class="line"><a name="l00663"></a><span class="lineno">  663</span>&#160;    </div>
<div class="line"><a name="l00665"></a><span class="lineno">  665</span>&#160;    <span class="keyword">struct </span>SparseMatrixBlock {</div>
<div class="line"><a name="l00666"></a><span class="lineno">  666</span>&#160;        EIGEN_MAKE_ALIGNED_OPERATOR_NEW</div>
<div class="line"><a name="l00667"></a><span class="lineno">  667</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> val;</div>
<div class="line"><a name="l00668"></a><span class="lineno">  668</span>&#160;        Eigen::VectorXi innerIdx, outerIdx;</div>
<div class="line"><a name="l00669"></a><span class="lineno">  669</span>&#160;    } uuT;</div>
<div class="line"><a name="l00670"></a><span class="lineno">  670</span>&#160; </div>
<div class="line"><a name="l00673"></a><span class="lineno">  673</span>&#160;    <span class="keywordtype">void</span> compute_uuT() {</div>
<div class="line"><a name="l00674"></a><span class="lineno">  674</span>&#160;        Eigen::MatrixXi pos=Eigen::MatrixXi::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, -1);</div>
<div class="line"><a name="l00675"></a><span class="lineno">  675</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00676"></a><span class="lineno">  676</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++)</div>
<div class="line"><a name="l00677"></a><span class="lineno">  677</span>&#160;            <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator it(<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>, i); it; ++it)</div>
<div class="line"><a name="l00678"></a><span class="lineno">  678</span>&#160;                <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator jt(<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>, i); jt; ++jt)</div>
<div class="line"><a name="l00679"></a><span class="lineno">  679</span>&#160;                    pos(it.row(), jt.row())=1;</div>
<div class="line"><a name="l00680"></a><span class="lineno">  680</span>&#160; </div>
<div class="line"><a name="l00681"></a><span class="lineno">  681</span>&#160;        uuT.outerIdx.resize(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>+1);</div>
<div class="line"><a name="l00682"></a><span class="lineno">  682</span>&#160;        uuT.innerIdx.resize(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00683"></a><span class="lineno">  683</span>&#160;        <span class="keywordtype">int</span> <a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>=0;</div>
<div class="line"><a name="l00684"></a><span class="lineno">  684</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) {</div>
<div class="line"><a name="l00685"></a><span class="lineno">  685</span>&#160;            uuT.outerIdx(j)=<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>;</div>
<div class="line"><a name="l00686"></a><span class="lineno">  686</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; i++)</div>
<div class="line"><a name="l00687"></a><span class="lineno">  687</span>&#160;                <span class="keywordflow">if</span> (pos(i, j)!=-1) {</div>
<div class="line"><a name="l00688"></a><span class="lineno">  688</span>&#160;                    uuT.innerIdx(<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>)=i;</div>
<div class="line"><a name="l00689"></a><span class="lineno">  689</span>&#160;                    pos(i, j)=<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>++;</div>
<div class="line"><a name="l00690"></a><span class="lineno">  690</span>&#160;                }</div>
<div class="line"><a name="l00691"></a><span class="lineno">  691</span>&#160;        }</div>
<div class="line"><a name="l00692"></a><span class="lineno">  692</span>&#160;        uuT.outerIdx(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>)=<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>;</div>
<div class="line"><a name="l00693"></a><span class="lineno">  693</span>&#160;        uuT.innerIdx.conservativeResize(<a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>);</div>
<div class="line"><a name="l00694"></a><span class="lineno">  694</span>&#160;        uuT.val=MatrixX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>*4, <a class="code" href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">nnz</a>*4);</div>
<div class="line"><a name="l00695"></a><span class="lineno">  695</span>&#160; </div>
<div class="line"><a name="l00696"></a><span class="lineno">  696</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00697"></a><span class="lineno">  697</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++)</div>
<div class="line"><a name="l00698"></a><span class="lineno">  698</span>&#160;            <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator it(<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>, i); it; ++it)</div>
<div class="line"><a name="l00699"></a><span class="lineno">  699</span>&#160;                <span class="keywordflow">for</span> (<span class="keyword">typename</span> SparseMatrix::InnerIterator jt(<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>, i); jt; ++jt)</div>
<div class="line"><a name="l00700"></a><span class="lineno">  700</span>&#160;                    <span class="keywordflow">if</span> (it.row()&gt;=jt.row()) {</div>
<div class="line"><a name="l00701"></a><span class="lineno">  701</span>&#160;                        <span class="keywordtype">double</span> _w=it.value()*jt.value();</div>
<div class="line"><a name="l00702"></a><span class="lineno">  702</span>&#160;                        <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> _uuT(4*<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>, 4);</div>
<div class="line"><a name="l00703"></a><span class="lineno">  703</span>&#160;                        <a class="code" href="class_dem_1_1_dem_bones.html#a0c30c98fb363f888345d86e5a11acd14">Vector4</a> _u;</div>
<div class="line"><a name="l00704"></a><span class="lineno">  704</span>&#160;                        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> s=0; s&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>; s++) {</div>
<div class="line"><a name="l00705"></a><span class="lineno">  705</span>&#160;                            _u=<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(s, i).homogeneous();</div>
<div class="line"><a name="l00706"></a><span class="lineno">  706</span>&#160;                            _uuT.blk4(s, 0)=_w*_u*_u.transpose();</div>
<div class="line"><a name="l00707"></a><span class="lineno">  707</span>&#160;                        }</div>
<div class="line"><a name="l00708"></a><span class="lineno">  708</span>&#160;                        <span class="keywordtype">int</span> p=pos(it.row(), jt.row())*4;</div>
<div class="line"><a name="l00709"></a><span class="lineno">  709</span>&#160;                        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> c=0; c&lt;4; c++)</div>
<div class="line"><a name="l00710"></a><span class="lineno">  710</span>&#160;                            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> r=0; r&lt;4*<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>; r++)</div>
<div class="line"><a name="l00711"></a><span class="lineno">  711</span>&#160;                                #pragma omp atomic</div>
<div class="line"><a name="l00712"></a><span class="lineno">  712</span>&#160;                                uuT.val(r, p+c)+=_uuT(r, c);</div>
<div class="line"><a name="l00713"></a><span class="lineno">  713</span>&#160;                    }</div>
<div class="line"><a name="l00714"></a><span class="lineno">  714</span>&#160; </div>
<div class="line"><a name="l00715"></a><span class="lineno">  715</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; i++)</div>
<div class="line"><a name="l00716"></a><span class="lineno">  716</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=i+1; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++)</div>
<div class="line"><a name="l00717"></a><span class="lineno">  717</span>&#160;                <span class="keywordflow">if</span> (pos(i, j)!=-1)</div>
<div class="line"><a name="l00718"></a><span class="lineno">  718</span>&#160;                    uuT.val.middleCols(pos(i, j)*4, 4)=uuT.val.middleCols(pos(j, i)*4, 4);</div>
<div class="line"><a name="l00719"></a><span class="lineno">  719</span>&#160;    }</div>
<div class="line"><a name="l00720"></a><span class="lineno">  720</span>&#160; </div>
<div class="line"><a name="l00721"></a><span class="lineno">  721</span>&#160; </div>
<div class="line"><a name="l00722"></a><span class="lineno">  722</span>&#160; </div>
<div class="line"><a name="l00724"></a><span class="lineno">  724</span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> mTm;</div>
<div class="line"><a name="l00725"></a><span class="lineno">  725</span>&#160; </div>
<div class="line"><a name="l00728"></a><span class="lineno">  728</span>&#160;    <span class="keywordtype">void</span> compute_mTm() {</div>
<div class="line"><a name="l00729"></a><span class="lineno">  729</span>&#160;        Eigen::MatrixXi idx(2, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>+1)/2);</div>
<div class="line"><a name="l00730"></a><span class="lineno">  730</span>&#160;        <span class="keywordtype">int</span> nPairs=0;</div>
<div class="line"><a name="l00731"></a><span class="lineno">  731</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; i++)</div>
<div class="line"><a name="l00732"></a><span class="lineno">  732</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=i; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) {</div>
<div class="line"><a name="l00733"></a><span class="lineno">  733</span>&#160;                idx(0, nPairs)=i;</div>
<div class="line"><a name="l00734"></a><span class="lineno">  734</span>&#160;                idx(1, nPairs)=j;</div>
<div class="line"><a name="l00735"></a><span class="lineno">  735</span>&#160;                nPairs++;</div>
<div class="line"><a name="l00736"></a><span class="lineno">  736</span>&#160;            }</div>
<div class="line"><a name="l00737"></a><span class="lineno">  737</span>&#160; </div>
<div class="line"><a name="l00738"></a><span class="lineno">  738</span>&#160;        mTm=MatrixX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>*4);</div>
<div class="line"><a name="l00739"></a><span class="lineno">  739</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00740"></a><span class="lineno">  740</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> p=0; p&lt;nPairs; p++) {</div>
<div class="line"><a name="l00741"></a><span class="lineno">  741</span>&#160;            <span class="keywordtype">int</span> i=idx(0, p);</div>
<div class="line"><a name="l00742"></a><span class="lineno">  742</span>&#160;            <span class="keywordtype">int</span> j=idx(1, p);</div>
<div class="line"><a name="l00743"></a><span class="lineno">  743</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++)</div>
<div class="line"><a name="l00744"></a><span class="lineno">  744</span>&#160;                mTm.blk4(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k)*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>+i, j)+=<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.blk4(k, i).template topRows&lt;3&gt;().transpose()*<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.blk4(k, j).template topRows&lt;3&gt;();</div>
<div class="line"><a name="l00745"></a><span class="lineno">  745</span>&#160;            <span class="keywordflow">if</span> (i!=j) <span class="keywordflow">for</span> (<span class="keywordtype">int</span> s=0; s&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>; s++) mTm.blk4(s*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>+j, i)=mTm.blk4(s*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>+i, j);</div>
<div class="line"><a name="l00746"></a><span class="lineno">  746</span>&#160;        }</div>
<div class="line"><a name="l00747"></a><span class="lineno">  747</span>&#160;    }</div>
<div class="line"><a name="l00748"></a><span class="lineno">  748</span>&#160; </div>
<div class="line"><a name="l00750"></a><span class="lineno">  750</span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> aTb;</div>
<div class="line"><a name="l00751"></a><span class="lineno">  751</span>&#160; </div>
<div class="line"><a name="l00754"></a><span class="lineno">  754</span>&#160;    <span class="keywordtype">void</span> compute_aTb() {</div>
<div class="line"><a name="l00755"></a><span class="lineno">  755</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00756"></a><span class="lineno">  756</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++)</div>
<div class="line"><a name="l00757"></a><span class="lineno">  757</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++)</div>
<div class="line"><a name="l00758"></a><span class="lineno">  758</span>&#160;                <span class="keywordflow">if</span> ((aTb(j, i)==0)&amp;&amp;(ws(j, i)&gt;<a class="code" href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">weightEps</a>))</div>
<div class="line"><a name="l00759"></a><span class="lineno">  759</span>&#160;                    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=0; k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>; k++)</div>
<div class="line"><a name="l00760"></a><span class="lineno">  760</span>&#160;                        aTb(j, i)+=<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, i).template cast&lt;_Scalar&gt;().dot(<a class="code" href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">m</a>.blk4(k, j).template topRows&lt;3&gt;()*<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(<a class="code" href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">subjectID</a>(k), i).homogeneous());</div>
<div class="line"><a name="l00761"></a><span class="lineno">  761</span>&#160;    }</div>
<div class="line"><a name="l00762"></a><span class="lineno">  762</span>&#160; </div>
<div class="line"><a name="l00764"></a><span class="lineno">  764</span>&#160;    _Scalar modelSize;</div>
<div class="line"><a name="l00765"></a><span class="lineno">  765</span>&#160;    </div>
<div class="line"><a name="l00767"></a><span class="lineno">  767</span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#aa8ea3a6b166652e954decbf135b3255b">SparseMatrix</a> laplacian;</div>
<div class="line"><a name="l00768"></a><span class="lineno">  768</span>&#160; </div>
<div class="line"><a name="l00770"></a><span class="lineno">  770</span>&#160;    Eigen::SparseLU&lt;SparseMatrix&gt; smoothSolver;</div>
<div class="line"><a name="l00771"></a><span class="lineno">  771</span>&#160; </div>
<div class="line"><a name="l00774"></a><span class="lineno">  774</span>&#160;    <span class="keywordtype">void</span> computeSmoothSolver() {</div>
<div class="line"><a name="l00775"></a><span class="lineno">  775</span>&#160;        <span class="keywordtype">int</span> nFV=(int)<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>.size();</div>
<div class="line"><a name="l00776"></a><span class="lineno">  776</span>&#160; </div>
<div class="line"><a name="l00777"></a><span class="lineno">  777</span>&#160;        _Scalar epsDis=0;</div>
<div class="line"><a name="l00778"></a><span class="lineno">  778</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> f=0; f&lt;nFV; f++) {</div>
<div class="line"><a name="l00779"></a><span class="lineno">  779</span>&#160;            <span class="keywordtype">int</span> nf=(int)<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>[f].size();</div>
<div class="line"><a name="l00780"></a><span class="lineno">  780</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> g=0; g&lt;nf; g++) {</div>
<div class="line"><a name="l00781"></a><span class="lineno">  781</span>&#160;                <span class="keywordtype">int</span> i=<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>[f][g];</div>
<div class="line"><a name="l00782"></a><span class="lineno">  782</span>&#160;                <span class="keywordtype">int</span> j=<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>[f][(g+1)%nf];</div>
<div class="line"><a name="l00783"></a><span class="lineno">  783</span>&#160;                epsDis+=(<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.col(i)-<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.col(j)).norm();</div>
<div class="line"><a name="l00784"></a><span class="lineno">  784</span>&#160;            }</div>
<div class="line"><a name="l00785"></a><span class="lineno">  785</span>&#160;        }</div>
<div class="line"><a name="l00786"></a><span class="lineno">  786</span>&#160;        epsDis=epsDis*<a class="code" href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">weightEps</a>/(_Scalar)<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>;</div>
<div class="line"><a name="l00787"></a><span class="lineno">  787</span>&#160; </div>
<div class="line"><a name="l00788"></a><span class="lineno">  788</span>&#160;        std::vector&lt;Triplet, Eigen::aligned_allocator&lt;Triplet&gt;&gt; triplet;</div>
<div class="line"><a name="l00789"></a><span class="lineno">  789</span>&#160;        <a class="code" href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">VectorX</a> d=VectorX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00790"></a><span class="lineno">  790</span>&#160;        std::vector&lt;std::set&lt;int&gt;&gt; isComputed(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00791"></a><span class="lineno">  791</span>&#160; </div>
<div class="line"><a name="l00792"></a><span class="lineno">  792</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00793"></a><span class="lineno">  793</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> f=0; f&lt;nFV; f++) {</div>
<div class="line"><a name="l00794"></a><span class="lineno">  794</span>&#160;            <span class="keywordtype">int</span> nf=(int)<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>[f].size();</div>
<div class="line"><a name="l00795"></a><span class="lineno">  795</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> g=0; g&lt;nf; g++) {</div>
<div class="line"><a name="l00796"></a><span class="lineno">  796</span>&#160;                <span class="keywordtype">int</span> i=<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>[f][g];</div>
<div class="line"><a name="l00797"></a><span class="lineno">  797</span>&#160;                <span class="keywordtype">int</span> j=<a class="code" href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">fv</a>[f][(g+1)%nf];</div>
<div class="line"><a name="l00798"></a><span class="lineno">  798</span>&#160; </div>
<div class="line"><a name="l00799"></a><span class="lineno">  799</span>&#160;                <span class="keywordtype">bool</span> needCompute=<span class="keyword">false</span>;</div>
<div class="line"><a name="l00800"></a><span class="lineno">  800</span>&#160;<span class="preprocessor">                #pragma omp critical </span></div>
<div class="line"><a name="l00801"></a><span class="lineno">  801</span>&#160;                <span class="keywordflow">if</span> (isComputed[i].find(j)==isComputed[i].end()) {</div>
<div class="line"><a name="l00802"></a><span class="lineno">  802</span>&#160;                    needCompute=<span class="keyword">true</span>;</div>
<div class="line"><a name="l00803"></a><span class="lineno">  803</span>&#160;                    isComputed[i].insert(j);</div>
<div class="line"><a name="l00804"></a><span class="lineno">  804</span>&#160;                    isComputed[j].insert(i);</div>
<div class="line"><a name="l00805"></a><span class="lineno">  805</span>&#160;                }</div>
<div class="line"><a name="l00806"></a><span class="lineno">  806</span>&#160; </div>
<div class="line"><a name="l00807"></a><span class="lineno">  807</span>&#160;                <span class="keywordflow">if</span> (needCompute) {</div>
<div class="line"><a name="l00808"></a><span class="lineno">  808</span>&#160;                    <span class="keywordtype">double</span> val=0;</div>
<div class="line"><a name="l00809"></a><span class="lineno">  809</span>&#160;                    <span class="keywordflow">for</span> (<span class="keywordtype">int</span> s=0; s&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>; s++) {</div>
<div class="line"><a name="l00810"></a><span class="lineno">  810</span>&#160;                        <span class="keywordtype">double</span> du=(<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(s, i)-<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(s, j)).norm();</div>
<div class="line"><a name="l00811"></a><span class="lineno">  811</span>&#160;                        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> k=<a class="code" href="class_dem_1_1_dem_bones.html#adf3d4f2aa7408bb6a9db540e86fd69c7">fStart</a>(s); k&lt;<a class="code" href="class_dem_1_1_dem_bones.html#adf3d4f2aa7408bb6a9db540e86fd69c7">fStart</a>(s+1); k++)</div>
<div class="line"><a name="l00812"></a><span class="lineno">  812</span>&#160;                            val+=pow((<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, i).template cast&lt;_Scalar&gt;()-<a class="code" href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">v</a>.vec3(k, j).template cast&lt;_Scalar&gt;()).norm()-du, 2);</div>
<div class="line"><a name="l00813"></a><span class="lineno">  813</span>&#160;                    }</div>
<div class="line"><a name="l00814"></a><span class="lineno">  814</span>&#160;                    val=1/(sqrt(val/<a class="code" href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">nF</a>)+epsDis);</div>
<div class="line"><a name="l00815"></a><span class="lineno">  815</span>&#160; </div>
<div class="line"><a name="l00816"></a><span class="lineno">  816</span>&#160;<span class="preprocessor">                    #pragma omp critical</span></div>
<div class="line"><a name="l00817"></a><span class="lineno">  817</span>&#160;                    triplet.push_back(<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(i, j, -val));</div>
<div class="line"><a name="l00818"></a><span class="lineno">  818</span>&#160;<span class="preprocessor">                    #pragma omp atomic</span></div>
<div class="line"><a name="l00819"></a><span class="lineno">  819</span>&#160;                    d(i)+=val;</div>
<div class="line"><a name="l00820"></a><span class="lineno">  820</span>&#160; </div>
<div class="line"><a name="l00821"></a><span class="lineno">  821</span>&#160;<span class="preprocessor">                    #pragma omp critical</span></div>
<div class="line"><a name="l00822"></a><span class="lineno">  822</span>&#160;                    triplet.push_back(<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(j, i, -val));</div>
<div class="line"><a name="l00823"></a><span class="lineno">  823</span>&#160;<span class="preprocessor">                    #pragma omp atomic</span></div>
<div class="line"><a name="l00824"></a><span class="lineno">  824</span>&#160;                    d(j)+=val;</div>
<div class="line"><a name="l00825"></a><span class="lineno">  825</span>&#160;                }</div>
<div class="line"><a name="l00826"></a><span class="lineno">  826</span>&#160;            }</div>
<div class="line"><a name="l00827"></a><span class="lineno">  827</span>&#160;        }</div>
<div class="line"><a name="l00828"></a><span class="lineno">  828</span>&#160; </div>
<div class="line"><a name="l00829"></a><span class="lineno">  829</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++)</div>
<div class="line"><a name="l00830"></a><span class="lineno">  830</span>&#160;            triplet.push_back(<a class="code" href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Triplet</a>(i, i, d(i)));</div>
<div class="line"><a name="l00831"></a><span class="lineno">  831</span>&#160; </div>
<div class="line"><a name="l00832"></a><span class="lineno">  832</span>&#160;        laplacian.resize(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>, <a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>);</div>
<div class="line"><a name="l00833"></a><span class="lineno">  833</span>&#160;        laplacian.setFromTriplets(triplet.begin(), triplet.end());</div>
<div class="line"><a name="l00834"></a><span class="lineno">  834</span>&#160; </div>
<div class="line"><a name="l00835"></a><span class="lineno">  835</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++)</div>
<div class="line"><a name="l00836"></a><span class="lineno">  836</span>&#160;            <span class="keywordflow">if</span> (d(i)!=0) laplacian.row(i)/=d(i);</div>
<div class="line"><a name="l00837"></a><span class="lineno">  837</span>&#160; </div>
<div class="line"><a name="l00838"></a><span class="lineno">  838</span>&#160;        laplacian=<a class="code" href="class_dem_1_1_dem_bones.html#a26babdede07e8e220002e420b9161aa9">weightsSmoothStep</a>*laplacian+<a class="code" href="class_dem_1_1_dem_bones.html#aa8ea3a6b166652e954decbf135b3255b">SparseMatrix</a>((VectorX::Ones(<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>)).asDiagonal());</div>
<div class="line"><a name="l00839"></a><span class="lineno">  839</span>&#160;        smoothSolver.compute(laplacian);</div>
<div class="line"><a name="l00840"></a><span class="lineno">  840</span>&#160;    }</div>
<div class="line"><a name="l00841"></a><span class="lineno">  841</span>&#160; </div>
<div class="line"><a name="l00843"></a><span class="lineno">  843</span>&#160;    <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a> ws;</div>
<div class="line"><a name="l00844"></a><span class="lineno">  844</span>&#160; </div>
<div class="line"><a name="l00847"></a><span class="lineno">  847</span>&#160;    <span class="keywordtype">void</span> compute_ws() {</div>
<div class="line"><a name="l00848"></a><span class="lineno">  848</span>&#160;        ws=<a class="code" href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">w</a>.transpose();</div>
<div class="line"><a name="l00849"></a><span class="lineno">  849</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00850"></a><span class="lineno">  850</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j=0; j&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j++) ws.col(j)=smoothSolver.solve(ws.col(j));</div>
<div class="line"><a name="l00851"></a><span class="lineno">  851</span>&#160;        ws.transposeInPlace();</div>
<div class="line"><a name="l00852"></a><span class="lineno">  852</span>&#160; </div>
<div class="line"><a name="l00853"></a><span class="lineno">  853</span>&#160;<span class="preprocessor">        #pragma omp parallel for</span></div>
<div class="line"><a name="l00854"></a><span class="lineno">  854</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> i=0; i&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">nV</a>; i++) {</div>
<div class="line"><a name="l00855"></a><span class="lineno">  855</span>&#160;            ws.col(i)=ws.col(i).cwiseMax(0.0);</div>
<div class="line"><a name="l00856"></a><span class="lineno">  856</span>&#160;            _Scalar si=ws.col(i).sum();</div>
<div class="line"><a name="l00857"></a><span class="lineno">  857</span>&#160;            <span class="keywordflow">if</span> (si&lt;_Scalar(0.1)) ws.col(i)=VectorX::Constant(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, _Scalar(1)/<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>); <span class="keywordflow">else</span> ws.col(i)/=si;</div>
<div class="line"><a name="l00858"></a><span class="lineno">  858</span>&#160;        }</div>
<div class="line"><a name="l00859"></a><span class="lineno">  859</span>&#160;    }</div>
<div class="line"><a name="l00860"></a><span class="lineno">  860</span>&#160; </div>
<div class="line"><a name="l00862"></a><span class="lineno">  862</span>&#160;    ConvexLS&lt;_Scalar&gt; wSolver;</div>
<div class="line"><a name="l00863"></a><span class="lineno">  863</span>&#160; </div>
<div class="line"><a name="l00868"></a><span class="lineno">  868</span>&#160;    <span class="keywordtype">void</span> compute_aTa(<span class="keywordtype">int</span> i, <a class="code" href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">MatrixX</a>&amp; aTa) {</div>
<div class="line"><a name="l00869"></a><span class="lineno">  869</span>&#160;        aTa=MatrixX::Zero(<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>, <a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>);</div>
<div class="line"><a name="l00870"></a><span class="lineno">  870</span>&#160;        <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j1=0; j1&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j1++)</div>
<div class="line"><a name="l00871"></a><span class="lineno">  871</span>&#160;            <span class="keywordflow">for</span> (<span class="keywordtype">int</span> j2=j1; j2&lt;<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>; j2++) {</div>
<div class="line"><a name="l00872"></a><span class="lineno">  872</span>&#160;                <span class="keywordflow">for</span> (<span class="keywordtype">int</span> s=0; s&lt;<a class="code" href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">nS</a>; s++) aTa(j1, j2)+=<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(s, i).homogeneous().dot(mTm.blk4(s*<a class="code" href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">nB</a>+j1, j2)*<a class="code" href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">u</a>.vec3(s, i).homogeneous());</div>
<div class="line"><a name="l00873"></a><span class="lineno">  873</span>&#160;                <span class="keywordflow">if</span> (j1!=j2) aTa(j2, j1)=aTa(j1, j2);</div>
<div class="line"><a name="l00874"></a><span class="lineno">  874</span>&#160;            }</div>
<div class="line"><a name="l00875"></a><span class="lineno">  875</span>&#160;    }</div>
<div class="line"><a name="l00876"></a><span class="lineno">  876</span>&#160;};</div>
<div class="line"><a name="l00877"></a><span class="lineno">  877</span>&#160;    </div>
<div class="line"><a name="l00878"></a><span class="lineno">  878</span>&#160;}</div>
<div class="line"><a name="l00879"></a><span class="lineno">  879</span>&#160; </div>
<div class="line"><a name="l00880"></a><span class="lineno">  880</span>&#160;<span class="preprocessor">#ifdef DEM_BONES_DEM_BONES_MAT_BLOCKS_UNDEFINED</span></div>
<div class="line"><a name="l00881"></a><span class="lineno">  881</span>&#160;<span class="preprocessor">#undef blk4</span></div>
<div class="line"><a name="l00882"></a><span class="lineno">  882</span>&#160;<span class="preprocessor">#undef rotMat</span></div>
<div class="line"><a name="l00883"></a><span class="lineno">  883</span>&#160;<span class="preprocessor">#undef transVec</span></div>
<div class="line"><a name="l00884"></a><span class="lineno">  884</span>&#160;<span class="preprocessor">#undef vec3</span></div>
<div class="line"><a name="l00885"></a><span class="lineno">  885</span>&#160;<span class="preprocessor">#undef DEM_BONES_MAT_BLOCKS</span></div>
<div class="line"><a name="l00886"></a><span class="lineno">  886</span>&#160;<span class="preprocessor">#endif</span></div>
<div class="line"><a name="l00887"></a><span class="lineno">  887</span>&#160; </div>
<div class="line"><a name="l00888"></a><span class="lineno">  888</span>&#160;<span class="preprocessor">#endif</span></div>
</div><!-- fragment --></div><!-- contents -->
</div><!-- doc-content -->
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a3462ee570b8b809ecff8739000783813"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a3462ee570b8b809ecff8739000783813">Dem::DemBones::m</a></div><div class="ttdeci">MatrixX m</div><div class="ttdoc">Bone transformations, size = [4*nF*4, 4*nB], m.blk4(k, j) is the 4*4 relative transformation matrix o...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00139">DemBones.h:139</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a050b398499f5daa9e42119b0774f25d0"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a050b398499f5daa9e42119b0774f25d0">Dem::DemBones::rmse</a></div><div class="ttdeci">_Scalar rmse()</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00356">DemBones.h:356</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a153f756d2284f955e35229f15753264c"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a153f756d2284f955e35229f15753264c">Dem::DemBones::lockW</a></div><div class="ttdeci">VectorX lockW</div><div class="ttdoc">Skinning weights lock control, size = nV, lockW(i) is the amount of input skinning weights will be ke...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00134">DemBones.h:134</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aba10dbe999f33dfb5456c0f3554700a3"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aba10dbe999f33dfb5456c0f3554700a3">Dem::DemBones::weightsSmooth</a></div><div class="ttdeci">_Scalar weightsSmooth</div><div class="ttdoc">[parameter] Weights smoothness soft constraint, default = 1e-4</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00097">DemBones.h:97</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a8008b672232b4d2f4081e00eea2ff062"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a8008b672232b4d2f4081e00eea2ff062">Dem::DemBones::cbTransformationsEnd</a></div><div class="ttdeci">virtual void cbTransformationsEnd()</div><div class="ttdoc">Callback function invoked after each bone transformations update.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00391">DemBones.h:391</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aed55acae84af67a7e958c830b821e747"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aed55acae84af67a7e958c830b821e747">Dem::DemBones::iterTransformations</a></div><div class="ttdeci">const int &amp; iterTransformations</div><div class="ttdoc">[zero indexed, read only] Current bone transformations update iteration number that can be used for c...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00154">DemBones.h:154</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a277dac86c8474a12eb489bedf98ddee2"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a277dac86c8474a12eb489bedf98ddee2">Dem::DemBones::nV</a></div><div class="ttdeci">int nV</div><div class="ttdoc">Number of vertices, typically indexed by i.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00114">DemBones.h:114</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_af1fcf96f35f6c4662fd95998225b9bda"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#af1fcf96f35f6c4662fd95998225b9bda">Dem::DemBones::init</a></div><div class="ttdeci">void init()</div><div class="ttdoc">Initialize missing skinning weights and/or bone transformations.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00186">DemBones.h:186</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a610c7725152df6465cd5f49311a3eeed"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a610c7725152df6465cd5f49311a3eeed">Dem::DemBones::cbWeightsIterEnd</a></div><div class="ttdeci">virtual bool cbWeightsIterEnd()</div><div class="ttdoc">Callback function invoked after each local weights update iteration, stop iteration if return true.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00401">DemBones.h:401</a></div></div>
<div class="ttc" id="anamespace_dem_html"><div class="ttname"><a href="namespace_dem.html">Dem</a></div><div class="ttdef"><b>Definition:</b> <a href="_convex_l_s_8h_source.html#l00015">ConvexLS.h:15</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a3b125934b440ebcebb3ba44acda41366"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a3b125934b440ebcebb3ba44acda41366">Dem::DemBones::transAffineNorm</a></div><div class="ttdeci">_Scalar transAffineNorm</div><div class="ttdoc">[parameter] p-norm for bone translations affinity soft constraint, default = 4.0</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00090">DemBones.h:90</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a0c30c98fb363f888345d86e5a11acd14"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a0c30c98fb363f888345d86e5a11acd14">Dem::DemBones::Vector4</a></div><div class="ttdeci">Eigen::Matrix&lt; _Scalar, 4, 1 &gt; Vector4</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00074">DemBones.h:74</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aae8731033731fbe34df518e15137724a"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aae8731033731fbe34df518e15137724a">Dem::DemBones::u</a></div><div class="ttdeci">MatrixX u</div><div class="ttdoc">Geometry at the rest poses, size = [3*nS, nV], u.col(i).segment(3*s, 3) is the rest pose of vertex i ...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00128">DemBones.h:128</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a673287121735ea64e0e4c244bacc432b"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a673287121735ea64e0e4c244bacc432b">Dem::DemBones::computeTranformations</a></div><div class="ttdeci">void computeTranformations()</div><div class="ttdoc">Update bone transformations by running nTransIters iterations with transAffine and transAffineNorm re...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00237">DemBones.h:237</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a53f9859cd63a8c51f6f50fdfa7251681"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a53f9859cd63a8c51f6f50fdfa7251681">Dem::DemBones::transAffine</a></div><div class="ttdeci">_Scalar transAffine</div><div class="ttdoc">[parameter] Translations affinity soft constraint, default = 10.0</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00088">DemBones.h:88</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ac4f7e7cb1b5cdfd737ddfc9e682de6f5"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ac4f7e7cb1b5cdfd737ddfc9e682de6f5">Dem::DemBones::fv</a></div><div class="ttdeci">std::vector&lt; std::vector&lt; int &gt; &gt; fv</div><div class="ttdoc">Mesh topology, size=[number of polygons], fv[p] is the vector of vertex indices of polygon p.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00148">DemBones.h:148</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aef1b9e1d0815f32c8f8374376c8d2c22"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aef1b9e1d0815f32c8f8374376c8d2c22">Dem::DemBones::iterWeights</a></div><div class="ttdeci">const int &amp; iterWeights</div><div class="ttdoc">[zero indexed, read only] Current weights update iteration number that can be used for callback funct...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00157">DemBones.h:157</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aa8ea3a6b166652e954decbf135b3255b"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aa8ea3a6b166652e954decbf135b3255b">Dem::DemBones::SparseMatrix</a></div><div class="ttdeci">Eigen::SparseMatrix&lt; _Scalar &gt; SparseMatrix</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00076">DemBones.h:76</a></div></div>
<div class="ttc" id="anamespace_dem_html_a5f2dfa58c8c80cb8b494cbdf8d82f105"><div class="ttname"><a href="namespace_dem.html#a5f2dfa58c8c80cb8b494cbdf8d82f105">Dem::indexing_vector</a></div><div class="ttdeci">Eigen::CwiseNullaryOp&lt; indexing_functor_vector&lt; ArgType, IndexType &gt;, typename indexing_functor_vector&lt; ArgType, IndexType &gt;::VectorType &gt; indexing_vector(const Eigen::MatrixBase&lt; ArgType &gt; &amp;arg, const IndexType &amp;indices)</div><div class="ttdef"><b>Definition:</b> <a href="_indexing_8h_source.html#l00115">Indexing.h:115</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a14f83e468b4a1eb6e2c6c951e7b5a4d3"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a14f83e468b4a1eb6e2c6c951e7b5a4d3">Dem::DemBones::cbInitSplitBegin</a></div><div class="ttdeci">virtual void cbInitSplitBegin()</div><div class="ttdoc">Callback function invoked before each spliting of bone clusters in initialization.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00374">DemBones.h:374</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ac99a3d810ab4fc1286e9677ce5c322b0"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ac99a3d810ab4fc1286e9677ce5c322b0">Dem::DemBones::subjectID</a></div><div class="ttdeci">Eigen::VectorXi subjectID</div><div class="ttdoc">Subject index of the frame, size = nF, subjectID(k)=s, where fStart(s) &lt;= k &lt; fStart(s+1)</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00125">DemBones.h:125</a></div></div>
<div class="ttc" id="a_mat_blocks_8h_html"><div class="ttname"><a href="_mat_blocks_8h.html">MatBlocks.h</a></div><div class="ttdoc">Defines some macros to access sub-blocks of packing transformation/position matrices for convenience.</div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_adf3d4f2aa7408bb6a9db540e86fd69c7"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#adf3d4f2aa7408bb6a9db540e86fd69c7">Dem::DemBones::fStart</a></div><div class="ttdeci">Eigen::VectorXi fStart</div><div class="ttdoc">Start frame indices, size = nS+1, fStart(s), fStart(s+1) are data frames for subject s.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00123">DemBones.h:123</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a5d5a50f2d4622a4823065cdca03ac9aa"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a5d5a50f2d4622a4823065cdca03ac9aa">Dem::DemBones::cbTranformationsBegin</a></div><div class="ttdeci">virtual void cbTranformationsBegin()</div><div class="ttdoc">Callback function invoked before each bone transformations update.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00389">DemBones.h:389</a></div></div>
<div class="ttc" id="a_convex_l_s_8h_html"><div class="ttname"><a href="_convex_l_s_8h.html">ConvexLS.h</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aa7f14b41174dc22e81b089d0c5a56950"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aa7f14b41174dc22e81b089d0c5a56950">Dem::DemBones::nIters</a></div><div class="ttdeci">int nIters</div><div class="ttdoc">[parameter] Number of global iterations, default = 30</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00080">DemBones.h:80</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ad6ded1630b2bac8d192e6134823a37f8"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ad6ded1630b2bac8d192e6134823a37f8">Dem::DemBones::nnz</a></div><div class="ttdeci">int nnz</div><div class="ttdoc">[parameter] Number of non-zero weights per vertex, default = 8</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00095">DemBones.h:95</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a26babdede07e8e220002e420b9161aa9"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a26babdede07e8e220002e420b9161aa9">Dem::DemBones::weightsSmoothStep</a></div><div class="ttdeci">_Scalar weightsSmoothStep</div><div class="ttdoc">[parameter] Step size for the weights smoothness soft constraint, default = 1.0</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00099">DemBones.h:99</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a913edb0e3273c6b18012c000a0764915"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a913edb0e3273c6b18012c000a0764915">Dem::DemBones::nS</a></div><div class="ttdeci">int nS</div><div class="ttdoc">Number of subjects, typically indexed by s.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00118">DemBones.h:118</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aee5cdd95063d2b11d8a461f0e2efa9f9"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aee5cdd95063d2b11d8a461f0e2efa9f9">Dem::DemBones::VectorX</a></div><div class="ttdeci">Eigen::Matrix&lt; _Scalar, Eigen::Dynamic, 1 &gt; VectorX</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00073">DemBones.h:73</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a9afbf2acccc4126e1e76c2a41f3ac649"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a9afbf2acccc4126e1e76c2a41f3ac649">Dem::DemBones::nTransIters</a></div><div class="ttdeci">int nTransIters</div><div class="ttdoc">[parameter] Number of bone transformations update iterations per global iteration,...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00086">DemBones.h:86</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aa9cf74a32230ec9c61096a58cdc555bf"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aa9cf74a32230ec9c61096a58cdc555bf">Dem::DemBones::nF</a></div><div class="ttdeci">int nF</div><div class="ttdoc">Number of total frames, typically indexed by k, nF = fStart(nS)</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00120">DemBones.h:120</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aa5bb61bfa1a60d355f4cce7f9cc60d5c"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aa5bb61bfa1a60d355f4cce7f9cc60d5c">Dem::DemBones::DemBones</a></div><div class="ttdeci">DemBones()</div><div class="ttdoc">Constructor and setting default parameters.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00105">DemBones.h:105</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a9dc51d56fe59c8c3e1e221b9819360f9"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a9dc51d56fe59c8c3e1e221b9819360f9">Dem::DemBones::cbWeightsEnd</a></div><div class="ttdeci">virtual void cbWeightsEnd()</div><div class="ttdoc">Callback function invoked after each skinning weights update.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00386">DemBones.h:386</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a06065f12b701c9bb556b7a79d14e984b"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a06065f12b701c9bb556b7a79d14e984b">Dem::DemBones::iter</a></div><div class="ttdeci">const int &amp; iter</div><div class="ttdoc">[zero indexed, read only] Current global iteration number that can be used for callback functions</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00151">DemBones.h:151</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ac4698b8e39c770e51aea7122689edcdf"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ac4698b8e39c770e51aea7122689edcdf">Dem::DemBones::cbWeightsBegin</a></div><div class="ttdeci">virtual void cbWeightsBegin()</div><div class="ttdoc">Callback function invoked before each skinning weights update.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00384">DemBones.h:384</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a9404399e2528175da18fd80d7d4fc4f9"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a9404399e2528175da18fd80d7d4fc4f9">Dem::DemBones::Vector3</a></div><div class="ttdeci">Eigen::Matrix&lt; _Scalar, 3, 1 &gt; Vector3</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00075">DemBones.h:75</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ac6f69fb5e7a91330f21f48ad9f019212"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ac6f69fb5e7a91330f21f48ad9f019212">Dem::DemBones::Triplet</a></div><div class="ttdeci">Eigen::Triplet&lt; _Scalar &gt; Triplet</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00077">DemBones.h:77</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a8fa642adc3a5610f0a391fa5f437a9b3"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a8fa642adc3a5610f0a391fa5f437a9b3">Dem::DemBones::w</a></div><div class="ttdeci">SparseMatrix w</div><div class="ttdoc">Skinning weights, size = [nB, nV], w.col(i) are the skinning weights of vertex i, w(j,...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00131">DemBones.h:131</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a9079e21634edb0dd2502bfc3a1eedbd0"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a9079e21634edb0dd2502bfc3a1eedbd0">Dem::DemBones::nInitIters</a></div><div class="ttdeci">int nInitIters</div><div class="ttdoc">[parameter] Number of clustering update iterations in the initalization, default = 10</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00083">DemBones.h:83</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_afab4c23ab2a15d11225934c9e4252881"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#afab4c23ab2a15d11225934c9e4252881">Dem::DemBones::weightEps</a></div><div class="ttdeci">_Scalar weightEps</div><div class="ttdoc">[parameter] Epsilon for weights solver, default = 1e-15</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00101">DemBones.h:101</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a5e4b99c4a233b477116204f9f1fae9cd"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a5e4b99c4a233b477116204f9f1fae9cd">Dem::DemBones::MatrixX</a></div><div class="ttdeci">Eigen::Matrix&lt; _Scalar, Eigen::Dynamic, Eigen::Dynamic &gt; MatrixX</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00070">DemBones.h:70</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aa0b51ddfe06be09fdf3c74a00f5dd867"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aa0b51ddfe06be09fdf3c74a00f5dd867">Dem::DemBones::nB</a></div><div class="ttdeci">int nB</div><div class="ttdoc">Number of bones, typically indexed by j.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00116">DemBones.h:116</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ad2a018d01cdb2a44f8ee964ddb4e8278"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ad2a018d01cdb2a44f8ee964ddb4e8278">Dem::DemBones::Matrix4</a></div><div class="ttdeci">Eigen::Matrix&lt; _Scalar, 4, 4 &gt; Matrix4</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00071">DemBones.h:71</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a276046cc5a2a13802762d65c17ec67fb"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a276046cc5a2a13802762d65c17ec67fb">Dem::DemBones::v</a></div><div class="ttdeci">Eigen::Matrix&lt; _AniMeshScalar, Eigen::Dynamic, Eigen::Dynamic &gt; v</div><div class="ttdoc">Animated mesh sequence, size = [3*nF, nV], v.col(i).segment(3*k, 3) is the position of vertex i at fr...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00145">DemBones.h:145</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a978a27bd67c19d78bd36242a29876212"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a978a27bd67c19d78bd36242a29876212">Dem::DemBones::cbIterEnd</a></div><div class="ttdeci">virtual bool cbIterEnd()</div><div class="ttdoc">Callback function invoked after each global iteration update, stop iteration if return true.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00381">DemBones.h:381</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aaaf326326d6094637dd0a995a778241f"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aaaf326326d6094637dd0a995a778241f">Dem::DemBones::cbTransformationsIterEnd</a></div><div class="ttdeci">virtual bool cbTransformationsIterEnd()</div><div class="ttdoc">Callback function invoked after each local bone transformations update iteration, stop iteration if r...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00396">DemBones.h:396</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aeb37e92f6b19ed95f4d2793c5f3dd300"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aeb37e92f6b19ed95f4d2793c5f3dd300">Dem::DemBones::cbInitSplitEnd</a></div><div class="ttdeci">virtual void cbInitSplitEnd()</div><div class="ttdoc">Callback function invoked after each spliting of bone clusters in initialization.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00376">DemBones.h:376</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html"><div class="ttname"><a href="class_dem_1_1_dem_bones.html">Dem::DemBones</a></div><div class="ttdoc">Smooth skinning decomposition with rigid bones and sparse, convex weights.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00066">DemBones.h:66</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a8391d148fee49d4c056b3b62bc1c7c08"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a8391d148fee49d4c056b3b62bc1c7c08">Dem::DemBones::lockM</a></div><div class="ttdeci">Eigen::VectorXi lockM</div><div class="ttdoc">Bone transformation lock control, size = nB, lockM(j) is the amount of input transformations will be ...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00142">DemBones.h:142</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ac8d2821c730539490ad6ce956a0b62a6"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ac8d2821c730539490ad6ce956a0b62a6">Dem::DemBones::compute</a></div><div class="ttdeci">void compute()</div><div class="ttdoc">Skinning decomposition by nIters iterations of alternative updating weights and bone transformations.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00344">DemBones.h:344</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a3101275f535d5c2435747b8ea5406f7d"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a3101275f535d5c2435747b8ea5406f7d">Dem::DemBones::computeWeights</a></div><div class="ttdeci">void computeWeights()</div><div class="ttdoc">Update skinning weights by running nWeightsIters iterations with weightsSmooth and weightsSmoothStep ...</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00276">DemBones.h:276</a></div></div>
<div class="ttc" id="anamespace_dem_html_a86c6e47fcd8031262219c42b65cf84ed"><div class="ttname"><a href="namespace_dem.html#a86c6e47fcd8031262219c42b65cf84ed">Dem::indexing_row_col</a></div><div class="ttdeci">Eigen::CwiseNullaryOp&lt; indexing_functor_row_col&lt; ArgType, RowIndexType, ColIndexType &gt;, typename indexing_functor_row_col&lt; ArgType, RowIndexType, ColIndexType &gt;::MatrixType &gt; indexing_row_col(const Eigen::MatrixBase&lt; ArgType &gt; &amp;arg, const RowIndexType &amp;row_indices, const ColIndexType &amp;col_indices)</div><div class="ttdef"><b>Definition:</b> <a href="_indexing_8h_source.html#l00047">Indexing.h:47</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_aed480cb1bda1b33f20171f7f847c1b06"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#aed480cb1bda1b33f20171f7f847c1b06">Dem::DemBones::Matrix3</a></div><div class="ttdeci">Eigen::Matrix&lt; _Scalar, 3, 3 &gt; Matrix3</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00072">DemBones.h:72</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a6391078901046e3acbd582dba6654fdd"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a6391078901046e3acbd582dba6654fdd">Dem::DemBones::cbWeightsIterBegin</a></div><div class="ttdeci">virtual void cbWeightsIterBegin()</div><div class="ttdoc">Callback function invoked before each local weights update iteration.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00399">DemBones.h:399</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a18bcd81e4238a4a298c4ee4ca78ae2fb"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a18bcd81e4238a4a298c4ee4ca78ae2fb">Dem::DemBones::clear</a></div><div class="ttdeci">void clear()</div><div class="ttdoc">Clear all data.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00161">DemBones.h:161</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_ae381fbe6610c9f4b169b16106dfbae94"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#ae381fbe6610c9f4b169b16106dfbae94">Dem::DemBones::nWeightsIters</a></div><div class="ttdeci">int nWeightsIters</div><div class="ttdoc">[parameter] Number of weights update iterations per global iteration, default = 3</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00093">DemBones.h:93</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_af43507c4ae91a4bf6550bc2afc25c023"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#af43507c4ae91a4bf6550bc2afc25c023">Dem::DemBones::cbIterBegin</a></div><div class="ttdeci">virtual void cbIterBegin()</div><div class="ttdoc">Callback function invoked before each global iteration update.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00379">DemBones.h:379</a></div></div>
<div class="ttc" id="aclass_dem_1_1_dem_bones_html_a938532a49b8ae162b866f9b59a75a7ab"><div class="ttname"><a href="class_dem_1_1_dem_bones.html#a938532a49b8ae162b866f9b59a75a7ab">Dem::DemBones::cbTransformationsIterBegin</a></div><div class="ttdeci">virtual void cbTransformationsIterBegin()</div><div class="ttdoc">Callback function invoked before each local bone transformations update iteration.</div><div class="ttdef"><b>Definition:</b> <a href="_dem_bones_8h_source.html#l00394">DemBones.h:394</a></div></div>
<!-- start footer part -->
<div id="nav-path" class="navpath"><!-- id is needed for treeview function! -->
  <ul>
    <li class="navelem"><a class="el" href="dir_d44c64559bbebec7f509842c48db8b23.html">include</a></li><li class="navelem"><a class="el" href="dir_4df80cdcb6e25b788688d708f51a1dfc.html">DemBones</a></li><li class="navelem"><a class="el" href="_dem_bones_8h.html">DemBones.h</a></li>
    <li class="footer">Generated by
    <a href="http://www.doxygen.org/index.html">
    <img class="footer" src="doxygen.png" alt="doxygen"/></a> 1.8.16 </li>
  </ul>
</div>
</body>
</html>
